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ABSTRACT 

We describe a general method for modeling gamma-ray burst prompt emission, and 
determine the range of magnetic field strength, electron energy, Lorentz factor of the 
source, and the distance of the source from the central explosion that is needed to 
account for the prompt 7-ray emission of a typical long duration burst. We find that 
for the burst to be produced via the synchrotron process unphysical conditions are 
required - the distance of the source from the center of the explosion (ii^) must be 
larger than ~ lO^^cm and the source Lorentz factor >10^; for such a high Lorentz factor 
the deceleration radius {Rd) is less than even if the number density of particles 
in the surrounding medium is as small as ~ 0.1 cm~^. The result, Ry > Rd, is in 
contradiction with the early x-ray and optical afterglow data that show that 7-rays 
precede the afterglow flux that is produced by a decelerating forward shock. This 
problem for the synchrotron process applies to all long-GRBs other than those that 
have the low energy spectrum precisely In order for the synchrotron process to 

be a viable mechanism for long-bursts, the energy of electrons radiating in the 7-ray 
band needs to be continuously replenished by some acceleration mechanism during much 
of the observed spike in GRB lightcurve - this is not possible if GRB prompt radiation 
is produced in shocks (at least the kind that has been usually considered for GRBs) 
where particles are accelerated at the shock front and not as they travel down-stream 
and emit 7-rays, but might work in some different scenarios such as magnetic outflows. 

The synchrotron-self-Compton (SSC) process fares much better. There is a large 
solution space for a typical GRB prompt emission to be produced via the SSC process. 
The prompt optical emission accompanying the burst is found to be very bright 14 
mag; for 2; ~ 2) in the SSC model, which exceeds the observed flux (or upper limit) for 
most GRBs. The prompt optical is predicted to be even brighter for the sub-class of 
bursts that have the spectrum fi, oc z^" with a ~ 1 below the peak of vf,^. Surprisingly, 
there are no SSC solutions for bursts that have a ~ 1/3; these bursts might require 
continuous or repeated acceleration of electrons or some physics beyond the simplified, 
although generic, SSC model considered in this work. Continuous acceleration of elec- 
trons can also significantly reduce the optical flux that would otherwise accompany 
7-rays in the SSC model. 



Subject headings: gamma rays: bursts, theory, method: analytical - radiation mecha- 
nisms: non-thermal 
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Introduction 



The last 10 years have seen a rapid advance in our understanding of gamma-ray bursts, due 
mainly to the study of GRB afterglows. We now know that at least some of the long duration 
GRBs (that last for more than about 5s) are produced iii the collapse of massive (young) stars 
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We also have good estimates of the total energy and beaming for these explosions as well as the 
property of the medium within about 1 pc of the explosion (Rhoads 1999; Sari, Piran & Halpern, 
1999; Frail et al 2001; Panaitescu & Kumar, 2001). However, our understanding of how the prompt 
7-ray emission is generated, and the mechanism for energy transport from the central engine (via 
magnetic field or kinetic energy of protons-neutrons and/or electron-positron pairs) remains highly 
uncertain. 

The goal of this paper is to provide a nearly model independent way of modeling the prompt 
7-ray emission with synchrotron or synchrotron-self-Compton (SSC) processes. We determine the 
basic properties of the 7-ray source from the data, and then determine how these can be interpreted 
in currently popular models such as the internal/external shock model. 

In the next section we provide the basic idea and details of the technique we use to model 
7-ray emission (the idea in its early form can be found in Kumar et al. 2006), and in §3 & §4 we 
describe the synchrotron and SSC results, respectively. 



2. Modeling 7-ray emission: basic idea and technical formalism 

The starting point for our modeling of the prompt 7-ray emission in GRBs is the assumption 
that the radiation is produced via the synchrotron or synchrotron-self-Compton processeqj in a 
source moving relativistically outward from the inner engine. Figure [T] provides a cartoon descrip- 
tion of our model. For a simple GRB light-curve (LC) consisting of a single peak we determine the 
average source properties corresponding to the time when the observed light curve peaks, and for 
a multi-peak GRB LC our calculation applies to individual pulses or spikes in the lightcurve. 

The source property can be uniquely described by the following set of 5 parameters: the 



^Mechanisms such as the inverse-Compton scattering of "photospheric" emission from a hot fireball (cf. Lazzati 
et al. 2000; Broderick, 2005) are not modeled by the approach we have adopted. And if it were to turn out that the 
GRB prompt emission is produced by such a mechanism then the work presented here is of little relevance. 
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magnetic field strength (B) in Gauss, the number of radiating particles {N) i.e., electrons and 
positrons, the optical depth of the source to Thomson scattering (r), the Lorentz factor of the 
source with respect to the rest frame of the GRB host galaxy (F), and the minimum electron 
energjl^ ji at the location where particles are accelerated (all the variables we use in this paper are 
defined in table 1 for easy reference). In addition, the particle distribution above ji is taken to be 
a power-law function: dn/d'y oc 7^^. Particles cool as a result of radiative losses and with time, or 
as they travel away from the acceleration site, and the distribution function becomes steeper than 
the index p at some energy where radiative losses become important. We calculate the modified 
distribution self-consistently as discussed below. We constrain this 5D parameter space with at 
least 4 observed quantities ~ the vfi, peak frequency i^^, the spectral index below v-y, the fiux at 
u^, the decay time of a single pulse in a GRB LC t^; p, the power law index, is constrained by the 
high energy spectral index, for v > u^. 




Fig. 1. — A schematic representation of our model. Assuming that radiation is synchrotron and inverse 
Compton, the 7-ray source properties can be described by five parameters {"/i,T, B, N,t) that determine 
the observed flux at one instance in time. We take this time to be the peak of a pulse in a GRB lightcurve. 
All of the calculations presented in this work apply to one single pulse in a typical GRB prompt lightcurve, 
as shown in the top left corner. 



^The electron energy is jinieC^, however, for convenience we suppress the factor rUeC? . 
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Table 1. Definition of variables 



7, minimum electron LF in source comoving frame 

r bulk; LF of the source 

B magnetic field strength, in Gauss, in source comoving frame 

T optical depth to Thomson scattering 

N number of radiating electrons (isotropic equivalent) 

a The spectral index below the peak of lyf^, i.e. oc v" 

'Ya LF of electrons emitting synchrotron at Ua 

7c LF of electrons emitting synchrotron at Uc 

Tsh LF of the shocked gas wrt the unshocked gas 

ly-y observed peak frequency of GRB i^fi, spectrum (i^-y^ = i^^/lQ^cV) 

Vi synchrotron injection frequency in observer frame {vi^ = eV) 

Vc synchrotron cooling frequency in observer frame [v^^ = t'c/lO^ eV) 

fa synchrotron self absorption frequcuicy in observer frame 

SSC self absorption frequency, below which the fl'' spectral index is +1 
external medium wind parameter, n = {A/mp)r^^; = A/5 x lO^^g cm 

d^^^ luminosity distance in units of 10^* cm 

E± kinetic energy of electrons and positrons (lab frame; isotropic equivalent) 

Eb energy in magnetic field (lab frame; isotropic equivalent) 

£■53 isotropic equivalent of outflow energy in units of 10^^ ergs 

fs/ke Eb/E± - ratio of magnetic to energy 

fa synchrotron prompt optical flux (in R band, at 2 eV) 

fx synchrotron prompt x-ray flux, at 1 keV 

/-y observed flux (in mJy) at u-y 

synchrotron flux at peak - min(i/i, Vc) 

no density of circum-burst medium 

rie comoving electron density in unshocked shell 

p power law index of electron energy distribution 

FLy distance from ccnter-of-explosion at which the radiation is produced 

Rd deceleration radius 

tj duration of one pulse in GRB light-curve (observer frame) 

ta the time available for electrons to cool before being re-accelerated 

Y Compton parameter 

z redshift 
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A relativistic moving source of finite angular size 9j (as seen by an observer at the center of 
explosion) can be treated as spherically symmetric as long as < 6j. The angular size determined 
from afterglow modeling suggests that 6j is larger than about two degrees for all bursts for which we 
have good data (Frail et al. 2001; Panaitescu &i Kumar, 2001) and a number of lines of argument 
and evidence suggests that T is greater than about 100 (cf. Piran 1992; Lithwick & Sari, 2001). 
Therefore, we can treat the source for prompt 7-ray emission as spherically symmetric, and the 
numerical values we quote in this paper are all isotropic equivalent quantities; for instance N is the 
total number of radiating particles in the source assuming the source to be spherically symmetric. 



2.1. Synchrotron and inverse- Compton radiations: basic equations 

The synchrotron injection frequency, fj, corresponding to electron minimum energy 7j, is 

* 27rmec(l + z) ' 

(eg. Rybicki & Lightman, 1979; Wijers Sz Galama, 1999), where q is electron charge, rUe the electron 
mass, c the speed of light, and z is the burst redshift. The synchrotron cooling frequency, fc, the 
characteristic frequency at which electrons cooling on a time scale ta (observer frame) radiate, is 

187rgmec(l + z) 

""""453^2(1+^)2 

where ctt is the Thomson scattering cross-section, and Y the Compton parameter. 

For most of the calculations in this work we assume that electrons are accelerated only once, 
and the time scale for acceleration is taken to be much less than the duration of a pulse in the GRB 
light-curve (tj). One time acceleration is, for instance, believed to apply to shocks where electrons 
arc accelerated at the shock front (by crossing the front back and forth multiple times) and not while 
they travel downstream; the picture is likely very different in magnetic reconnection/dissipation. To 
capture some of the effects of multiple-times particle acceleration in time period of a pulse duration 
in GRB LC we introduce a time scale, ta, which is the average time in between two successive 
episodes of particle acceleration or the time available for electrons to cool in between acceleration; 
for one shot acceleration ta = tj, and in the opposite limit of continuous acceleration = 
when the rate of energy gain is balanced by radiative loss rate. The electron distribution function 
resulting from the acceleration process is taken to be dn/d'je oc -je ^ ■ The distribution function in 
the source as a whole is different due to the radiative cooling of electrons with time. The electron 
distribution function averaged over the source is described by two characteristic energies viz. 7^ 
and 7c; 7cmeC^ is the energy of electrons that cool on time scale ta- The electron distribution for 
7e > max(7j,7c) is proportional to , and the distribution between jc and ji (for 7c < 7i) is 

proportional to 7^^. Electrons cool via synchrotron and inverse-Compton losses. The rate of loss of 
energy is affected by the synchrotron self-absorption frequency Va — electrons with characteristic 



(1) 



(2) 



- 6 - 



synchrotron frequency below Va lose energy only via the inverse-Compt qn scattering process. W e 
calculate 7c and Ua by solving a coupled set of equations as described in lMcMahon et alj (|2006l ). 



The synchrotron flux at the peak of the fi, spectrum, at min [vi, f^], is given by 

VSq^BNF 



(3) 



where is the luminosity distance to the source. The effect of synchrotron self absorption is 
not included in the above expression for /j^^, and therefore the observed flux, in general, would 
be different from f^^. The flux at other frequencies are calculated as described in Sari, Piran and 
Narayan (1998). 

The inverse-Compton flux (in observer frame) at frequency u, f^'^{u), is calculated using the 
following equation (cf. Rybicki & Lightman, 1979) 



71 7e«7e 



(4) 



where f{vs) is the synchrotron flux per unit frequency in the observer frame, 5r is the radial extent 
of the source (comoving frame) which is related to the optical depth r (one of the five parameters 
we use to characterize the source), the function F{x) is 



F{x) = 2x\\ix X + \ — 2x^ , for < X < 1, 



(5) 



and 7i is the minimum LF for electron distribution. We include the Klein-Nishina correction to 
the above expression when Vgle/'^ > nieC? ■ 



The expression for the Compton Y parameter is 



Y = (TT 



dr' / d7e7e^^ = \r 



p — 1 
p-2 



if X \ 



^p{p-l) ^ 



7i ~ 7c P > 2 



(6) 



where the r'-integral is over the comoving radial width of the source. For ease of future use we 
rewrite the above expression for y as y = 4r7?^/3 where 



p-2 



X < 



(z^c/i^i)^ 7c«;7i P>2 

{uc/ly^f~''{3-p^^ ji-^^c 2<p<3 

{p-sy^ ji<s^^c p>3 



(7) 



^p{p- 1) 



-1 



7c P > 2. 



In our 5D parameter space search, we limit the Compton y-parameter for a synchrotron solution 
to be less than 10 and for an inverse-Compton solution y<100. The rationale for the constraint on 
the y-parameter is that we want an efficiency of >10% in the 7-ray energy band of ~ 10-400 keV; 
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observations suggest this efficiency for a typical long duration burst from a comparison of energy in 
the 7-ray radiation and the kinetic energy of the ejecta determined from the afterglow observations 
(Panaitescu & Kumar, 2002). 

We can calculate the distance of the source from the center of explosion with 2 of the 5 
parameters, A'" and r: 



We consider a 7-ray source of a finite lifetime, at a distance Fiy from the central explosion, 
that is responsible for generating one pulse in the observed GRB prompt lightcurve. Electrons in 
the source are heated during some time interval, set by the central engine variability/activity time, 
and subsequently the source undergoes adiabatic expansion. The width of an observed GRB pulse 
is determined by a number of different factors - the central engine variability time, the adiabatic 
expansion and cooling time, and the curvature time - which are described below. 

1. Central engine variability time: sets the observed GRB pulse width if the variability time 

is larger than ~ Ry/2cT'^ and if the distance where ^-ray photons are generated does not 
increase with time; in this case GRB pulse duration is independent of R^. However, when 
the source is turned off' the 'y-ray flux would decline on a timescale of Rj/2cT'^ , the adiabatic 
expansion time scale; for all the calculations in this work we use the light-curve decline time 
although for simplicity we continue to refer to it as pulse width. 

2. Curvature timescale: is the time interval between arrival of photons with angular separa- 
tion of F"^ as seen by an observer at the center of explosion. R is the minimum time scale 
for a 7-raj/ pulse width, as long as the outflow from the GRB has an angular size larger than 

, and is equal to R^/{cF'^). 

3. Adiabatic expansion time scale — this is the timescale for electrons/protons to cool because 
of expansion of the source. As the distance of the source from the center doubles its volume 
increases by a factor >4, and electron/proton energy drops by a factor 2. This timescale in 
the observer frame is ^ R^/{cF'^). 

We note that if electrons are heated by coupling with protons, then the time scale for electrons 
to cool down can be larger than R^/{cF'^); cooling time in this case can be as large as 
R^/ (cF^) X (proton energy/electron energy) provided that protons transfer energy to electrons 
for this time duration, and the energy transfer rate balances the loss of energy for electrons 
to adiabatic and radiative coolings. However, the coupling between p'^s and e^s is unlikely to 
increase the pulse width by a large factor (>5) unless the energy in electrons is much smaller 




(8) 



2.2. Relation between R^ and pulse-width 
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than in ]?+s, but in that case the efficiency of 7-ray generation would be small which is not 
supported by observations. 

4. Radiative cooling timescale. For Vc > ui the radiative cooling time scale is larger than the 
adiabatic time scale and in that case GRB pulse width is equal to the adiabatic or curvature 
time. For < Vi electrons cool on a smaller time scale, and once electron heating stops, 
the lightcurve falls off on the curvature timescale. We note that the observed pulse duration 
cannot be larger than i?-y/(2cr^) which corresponds to an elapsed time of R-^jc in the center 
of explosion frame, and during this time the source has moved to iiy from the center. 

Thus we see that the observed decay time for a pulse in GRB LC, produced by a relativistic source, 
is 

R^{\ + z) 
~ 2cr2 ' 

when the GRB redshift is z. 



2.3. Energy etc. 

The total energy in the source consists of the kinetic energy of electrons and positrons {Ej^ 
and the magnetic field [Eb)'- 

E± = N{p - l)^ime(?V/{p - 2), Eb = R'^B^ /Q. (10) 

Note that {p— l)^imcc'^ / (p — 2) is the average energy per electron/positron in the source comoving 
frame at the acceleration site, and in the calculation of Eb we took the comoving radial thickness 
of the source to be Ry /V which is roughly what one expects for a causally connected source where 
the signal speed is close to the speed of light. 

We do not make any assumptions regarding the energy in protons since protons do not con- 
tribute to the observed 7-ray radiation. This has the effect that the parameter space we determine 
is larger than it would be if protons carried a substantial amount of energy since the energy avail- 
able to would be smaller than the upper limit of 10^^ erg (isotropic equivalent) we impose in 
our search for solutions in the 5D parameter space. 

For a 7-ray source that arises from shock heated gas, the minimum electron energy behind 
the shock front, 7^ (one of the five parameters we use), can be related to the Lorcntz factor of the 
shocked gas wrt the unshocked gas, Tg^. The minimum Vg^ needed to produce 7^ is 



me{p - 1) 
_mp{p-2) 



27i, (11) 



where rup (me) is proton (electron) mass. The factor 2 in the above expression is for the case where 
there is an energy equipartition between electrons and protons and there are no pairs in the 
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plasma; will be larger if there are pairs or if electrons have less energy than protons O 

2.4. The basic technique for finding source properties 

We determine the properties of the 7-ray source for a GRB by finding the region in the 5-D 
parameter space {'yi,T,B,N,T) that satisfies the following set of observational constraints: 1) 
the frequency at the peak of the vfu spectrum (2^7); 2) the peak flux at v^; 3) the spectral index 
above z^^ - which constrains electron index p - and the index below v^; 4) the burst duration - 
for a GRB with a single pulse in the LC - or the duration of an individual pulse (t^) for GRBs 
with complicated LC; 5) optical and x-ray prompt flux or limit if available. The flux at a given 
observer time reflects the property of the source averaged over equal-arrival-time volume, therefore, 
the observed peak flux depends on the evolution of the source and this introduces uncertainty in 
the flux calculation by a factor of about two. For this reason we only require the theoretical flux 
to match the observed value to within a factor of ~ 2. 

We now use this technique to find the 5D solution space and source property for GRBs produced 
via synchrotron (§3) and SSC (§4). 

3. Synchrotron solutions 

We consider in this section the parameter space of solutions when the observed 7-rays are 
produced via the synchrotron procescl. First, we determine approximate solutions by analytically 
solving a system of equations for our 5 parameters (7^, F, B, N, r) for the generic synchrotron case. 
The solutions for each parameter are expressed in terms of the Compton Y parameter and three 
observed quantities: the frequency where I'fi, peaks, the 7-ray flux at this frequency (Z^; in 
mJy), and the duration of a pulse in GRB LC (t^); y is a convenient and useful parameter because 
its value is expected to lie in a limited range, e.g. Y<1 for the synchrotron solutions & l<y<10 for 
the SSC process. Having the general synchrotron solutions in hand, we then apply the analytical 
results to the low energy spectral index cases of a = 1/3, a = —1/2 and a = —{p — l)/2, compare 



^If there are "cold" protons and electrons in the shocked gas, i.e. only a fraction of particles in the shocked gas are 
accelerated, and electrons have more than mpC^Fgh energy, we would in that case overestimate r^h- However, in this 
work we do not consider that there is a cold component to the 7-ray source since such a component would not radiate 
and affect observations and the solutions in the 5-D parameter space. Therefore, such a cold component, if present, 
would have little effect on all of the major results in this work; the only quantity affected by the cold component is 
the value for Tsh which is a peripheral quantity and not part of the central flow of the logic in this paper. 

* It has been suggested that another radiation process, called jitter, might be responsible for 7-ray generation for 
those bursts that have low energy spectrum f^ocu (Medvedev, 2000). We show in appendix B that whenever jitter 
radiation dominates the observed flux to produce a. fu ot v spectrum the Compton-Y parameter is extremely large - 
F>10'' - and most of the energy of the explosion comes out in ~ 100 GeV SSC photons. 
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the analytical and numerical results, and draw conclusions as to the process by which 7-rays are 
generated in GRBs; the spectral index a is defined by /jy oc for u < v^. 

The 5 equations that we solve are those for the observer frame synchrotron injection frequency 
Ui ([H), the cooling frequency Uc ^1 the pulse duration (i^) ([8]&[9j), synchrotron flux f^^ in mJy at 
Vp = min(fj, Uc) ([3]), and the Compton Y parameter ([6]): 

u,, = 1.1 X W-'^B^fr{l + z)-' (12) 

i/c5 = 6.6 X 10^(1 + z)B^^r^^t^'^{l + y)~2 (13) 

r pa 1.5 X 10^A^55(1 + zfr^H-^ (14) 

= IWBTN^^dll (1 + z) mJy (15) 

where = Ui/lO^eV, = z^c/lO^eV, and iYgg = iV/10^^ 

To solve Equations [T2]fT5] & [H we first eliminate A^55 from equation (fT5j) using equation (fH|) , 
then eliminate r using equation ([6]) to find 

srSr' « 1-9 X ioV.,dLs(i + ^Y^^y^i- (16) 

Next, combining ([TB]) & (fT2|) we get 

7-4r4 ^ 2.1 X 10- V,rV,^4^^y-iet-^. (17) 

Multiplying the square root of equations (jl2p and (jlSj) together, we have 

B^r^ « 8.5 X 10-5t-iz.,7i/,7(l + y)-^ (18) 

We can eliminate 7^ from equations (fT7|) and ([TS]) by dividing equation (fTTl) by ([T8|) to the fourth 
power: 

T'B-^ « 4.0 X + Yfitt-^ (19) 

And finally, if we multiply equation (jl9p by the fourth power of (jl6p and divide by the square of 
(|17p we find the solution for F to be 

r « io3z/|4/|t;^i|i^-^(i + 5^)^^428(1 + ^)^- (20) 

Using r, we can solve for "fi, B, and r: 

ji^4.7xW''u^'ulj~-^H'^tiY^{l+Y)k-^dll{l + z)-* (21) 

S « 4.0i/,r^i.e7/.;^4*a^^^(l + ^)^^r^dZi(l + ^)^ Gauss (22) 
T^3.3xlO-''>u~K-^fltyhahl{l+Y)-H-ldlJl + z)'-2. (23) 



- 11 - 



Equations (j20p -(|23 p provide approximate solutions for {'yi,T,B,N,T) when the synchrotron pro- 
cess produces the observed 7-ray radiation; more accurate solutions for these parameters are ob- 
tained by numerical calculations and the results are shown in Figures [2H8l These general solutions 
can be used to investigate different cases of low energy spectral indices (a) by adopting appropriate 
values for and v^- The full dependences on these two frequencies are not completely shown 
here - each case of a has a different functional dependence on ^, and ^ is a function of Vi and Vc- 

Note that f^^ is not the observed flux at z^^, the peak of 7-ray spectrum, but is the flux at 
min(t'j, fc) = Vp, and the effect of synchrotron-self-absorption, if any, at Vp has been ignored. Since 
the dependence of the parameters F, etc. on f^^ is very weak (eqs. [20H23]) . we do not worry 
about the difference between 7^,^ & at this point, even though fi^^ can be much greater than 
(the flux at u^); is the peak of vfy - for p < 3, = max[^'j, Vc\ & for p > 3 = min[fj, Uc\). 

Using the parameter solutions, we can derive the distance of the 7-ray source from the center of 
explosion (-R-y), and the energy in the magnetic field and electrons. The radius i?-y = 2cF^t-y(l-|-z)~^ 
is found to be 

^ 6.0 X lQ^^ul/}Jl^tlY~i{l + Y)hildlJl + z)~'2 cm (24) 

and should be compared to the deceleration radius {Rd) of the GRB outflow in both a homoge- 
neous external medium with particle number density ng, and a wind external medium where the 
particle number density is given by {A/mp)r~'^ (these are two special cases of a power law density 
stratification - the density varying as r~* - corresponding to s = & s = 2) 

(25) 



Ra 



1.2 X 10i'^£;|3no ■^F2 ' cm s = 



1.8 X lO^^^.v^^^^F^^ cm 



where £^53 is the isotropic equivalent energy in GRB-ejecta in units of 10^^ ergs, F2 = F/lOO, and 
Aj, = A/{f) X lO^^g cm~^). Substituting in the solution for F, we find R^ to be 



R. 



1111 



(26) 



2.6 X l{)^^El^n^-'u^^y^^^^ Uj^f^ta^'Yiil + Y)-iC^^d^l{l + z)-i cm s = 

l.SxlQ^^E5^A-^u~}u~}u}^tf^Yl{l + Y)-'2^-ldll{l + z)-'2 cm s = 2 

The magnetic and energies, given by equation (jlOp . are found to be 

Eb ~ 5.8 X 10^0 i^iJu,hY-^Cdl,,il + z)-' ergs (27) 

and 

E± « 8.5 X 10^° 4''cJuM'^ + ^)4.s(l + ^)"' ergs. (28) 

Since the dependence of the above two quantities on f^^ is linear, we should replace f^^ with 
the flux observed at i/^. This will be done in the following sections, since the expression for 
depends on a. 
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We can relate the solution-subspace we find to parameters for the shock model for GRBs; if 
electrons are accelerated in a relativistic shock then the LF of shock front (Tsh) wrt to the unshocked 
material is related to ji (one of the 5 parameters) and is given by equation (jlip 

where is the ratio of energy in electrons and the total thermal energy in the 7-ray source. 

We now apply the results obtained in this section to each possible synchrotron low energy 
spectral index a. 



where 



3.1. Synchrotron solutions when the low energy spectrum is ly 2 

We use equation d?]) to eliminate ^ from the analytical solutions given by equations [20VI23I for 

the 7i <C 7c & 2 < p < 3 case, and substitute = & /j^^ = f^{u^^/ui^) 2 , to find that 
synchrotron solutions for a = — ^^2^^ 

19-3p JL _3 3p-9 1 1 1 3 _3_ 

r«io3z.^3«2 /«i^«z../^ t|y-T6(i + y)4(i + z)4d|^^A« (so) 

7i « 4.7 X 10^1.^32 16^8^.^32 tjY-e{l + y)4(l + z)4dJ^A,;« (31) 
B ^ 4.0u^ f~^t^u^ta^YTE{l + y)-t(l + z)UllA~,^ Gauss (32) 

15p-47 1 _1 29-15p _1 y , , 1 __7 

r«3.3x 10-i%^« f^h'^^,'' ta~'Y-^{l + Y)-2{l + zr2dl^^A^;' (33) 

and we should emphasize that is not the burst duration - it is the width of a single spike in the 
GRB prompt-lightcurve. 

For a typical long duration GRB with = ImJy, = lOOkeV, = 0.1s, z = 1, d^jg = 2, 
and ta ^ - henceforth we will refer to a GRB with these observed parameters as GRB-\] - the 
5-parameters of the 7-ray source {^i,T,B,N,T) are obtained from equations I30fl33] and are given 
by 

r> 3.2 X io3y-^(i + y)^ (35) 
7i< 6.0 X io3yiii(i + y)3 (36) 

5<16yi^(l + y)-iGauss (37) 
r> 4.7 X 10-^yi(l + y)-i (38) 

In deriving these inequalities we took p = 2.5, z^c = i/^ = 100 keV, and < 0.1. 
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The dependence of T, 7^, i? & r on y is weak, so the coefficients in above expressions are 
reasonable estimates for the 7-ray source basic physical parameters for GRB-\\. We see that the 
7-ray source LF, F, is required to be rather large - r>3 x 10^ - if the radiation is to be produced 
via the synchrotron process. This large F is not consistent with afterglow modeling, which gives 
a value of a few hundred or less (Panaitescu & Kumar, 2002). Furthermore, as shown below, the 
distance of 7-ray source from the center of explosion turns out to be larger than the deceleration 
radius for this large F value, unless no is very small. This suggests that the synchrotron solution 
is internally inconsistent; after the deceleration radius F is a function of A^, 7i and no and is no 
longer an independent parameter as considered in these derivations. The possibility that > Rd 
is also ruled out by early optical afterglow data - eg. GRBs 050801, 050820A, 060124, 060418, 
060607A, 060614, 060714 - that show that 7-rays precede a rising afterglow flux that is produced 
by a decelerating forward shock. Moreover, if then in this case of a decelerating source we 

should see an increasing GRB pulse duration with time, which is not observed. 

The distance of the 7-ray source from the center of explosion, R^ w 2c^F^/(l-|-z), is calculated 
using eq. (f30l) . and is given by 

19-3p 3 1 3p-9 13 1 1 3 3 

R^ « 6.0 X 10^^ u^,^^ /l^^is" tlY—^{l + y)2 (1 + zY^dl^^Al^ cm (39) 

or i?^ ~ 3 X 10^*^ y-i(l + y)5 cm for GRB\ 

We now compare these analytical estimates to the numerically computed solution space for 
synchrotron radiation. A numerical search of the allowed region of the 5-D parameter space that 
satisfies the observational constraints (i/^, & t^; the same constraints that we used in the deriva- 
tion of analytical expressions), confirms that for synchrotron solutions i?'y^l0 cm, F>10^ and 
10<7i<10^ (see fig. [2]). We have considered a wide range of values of peak frequency iy^)-, 7-ray 
flux at the peak, and pulse duration, to see if we can find some viable synchrotron solutions for any 
GRBs with a = — (p — l)/2. These solutions are shown in Figure [2j We find that by decreasing 
any of the observable parameters R-y decreases, but the dependence is weak in agreement with the 
scaling given in equation ([39|) . Furthermore, a decrease in reduces F as expected from equation 
(f30j) . but even for = 10 ms, F is stiU >l{f'. 

We next calculate the deceleration radius and compare it with R^ to ensure R^ < Rd for self 
consistent solutions. The deceleration radius for GRB-ejecta is calculated using eq. ([26]) and is 
given by 

1 _1 3p-19 _1 1 3-p _1 111 _i _1 

2.6 X 10^'' E^^Hq 3jy^g48 ta''Y8{l + y)-6 (1 + zyed^^^A^p' cm s = 

Rd-{ (40) 

3p-19 _ 3 3 9-3p _1_ 3 11 _3 _3 

[1.8 X 1013 E^sA-^u^,'' H^u,^'' ta ^Y-s{l + y)-2 (1 + z)-2ciJ^A,/ cm s = 2 
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- - — 400 keV 
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Fig. 2. — Results of numerical calculation for the allowed synchrotron solution space when the 
spectrum below the peak of vf^, at i/^, is: /jy oc for v < v^. A point in the 5-D parameter 

space {jijT, B, N,t) is considered an allowed solution for the observed GRB parameters {v.y, f-y, 
t^, a) provided that is within a factor 2 of the observed value, the pulse duration (t^) &; flux 
at (fy) are within a factor 1.5 &: 3 of the observed value respectively; the larger tolerance on 
flux is due to larger error in flux calculation. The x-axis shows the distance of the 7-ray source 
from the center of the explosion. The top left panel is ji - the minimum LF of electrons in source 
comoving frame at the site where they are accelerated (electron distribution function for > 7, is: 
drie/d'^e oc 7,r^'^ P = 2.5). The top right panel shows the bulk LF of the source, the bottom left 
panel shows the comoving magnetic field in Gauss, and the bottom right panel shows the ratio of 
energy in the magnetic field and electrons. For all of the numerical calculations we took the burst 
redshift z = \. Legend shows several different cases of GRBs corresponding to different observed 
values for v^, f^, and t^. Only one observational parameter - that noted in the legend - is changed 
at a time, all the remaining parameters are left unchanged; the base value for the parameters is the 
same as we took for GRB-\\, i.e. = 100 keV, = ImJy, t-y = 0.1s, and ta = t^. For instance, 
for the 20 keV case, denoted by the solid black line, i^^ = 20keV, and f-y &: t-y are same as for 
GRB-\\ i.e. 1 mJy and 0.1 s respectively. 
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and the ratio of and is: 



Rj 
R^ 



2.3i5;53^4z.^'-'/|i/"t|y-^(i + y)l(i + z^Ul^^aI s = 

19 -3p 3 _1 3p-9 3 3 3 

3.3xW^E^,'A,.u,,' f-H,-^v^^^' taY-^{l + Y)dlAl 



(41) 



Substituting in the observable parameters for GRB-\\ into the above equation and solving for no 
& At, such that R^/R^ < 1, we find 

no < 0.057£;53t-2yi(l + y)-2 cm-3 s = (42) 

A, < 5.5 X IQ-^E^-it-^Yi (1 + s = 2 (43) 

Note that no &: A,,, must be very small to ensure that R-y < Rd, especially for F < 1 expected of 
synchrotron solutions. Figure ([3]) shows the results of numerical calculations which confirms these 
analytical estimates. Moreover, if we want R.y/Rd'^0.5, in order to have a clear separation between 
internal and external shocks, then no^lO^^cm^^. Therefore, self-consistent synchrotron solutions 
with Ry^Rd require very low dens ity for the circumstellar me dium compared with no ~ 1 cm^'^ 



obtained from afterglow modeling (jPanaitescu &: Kumaij 120011 ). The limit on no can be increased 



by decreasing ta (see eq. H2|) . Numerical result for the upper limit on no when ta = t^/WO is 
shown in fig. [3l It confirms the analytical result that no ~ 1 cm~^ can give i?-y < Rd provided 
that ta <^ t^. It should be noted that for systems involving shock heating of particles we expect 
ta ~ ty because electrons are accelerated at the shock front and there is no subsequent acceleration 
as particles travel downstream; in magnetic reconnections or dissipation it is natural to expect 

We now estimate the ratio of energy in and magnetic field to find out if it is much less than 
unity or not when ta < t^ (a small value for E±/Eb results in low efficiency for 7-ray generation). 
The ratio E±/Eb can be calculated using eqs. (p7|) & ([28]) and is given by 



E± _ ^ f T^To 



+ Y), (44) 



for 2 < p < 3 (numerical calculations take p = 2.5). For the solution space corresponding to 
a = —(jp — l)/2, O.l^z^^/t'j^lO^ and so E± > Eb even when ta/t-y ~ 10"^. Therefore, small ta/t-y 
solutions are fine from the point of radiative efficiency; the above equation needs to be modified, 
when ta <C i-Y, to include the total energy input in electrons during a GRB pulse width of ty, which 
will further improve the radiative efficiency when ta/t-y is very small. 

The reason that these synchrotron solutions have large Ry is not hard to understand. It 
requires a certain minimum number of electrons to produce the observed flux of /-^ ~ 1 mJy at 
Uy ~ 100 keV: N ~ 10^^/(-Br) - see eq. [T5j And in order to keep the Compton-Y parameter, 
Y ~ T'ji'jc, less ~ 10 — otherwise most of the energy will come out in IC-scattered photons at ^ 
1 MeV — we must have large Ry for the source. The solution offered by ta <C ty is also easy to 
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[cm] [cm] 

Fig. 3. — Left panel: the upper limit to the ISM density (no) for synchrotron solutions with a = — (p — l)/2 
for a burst with v-y ~ 100 keV, f-y = I mJy, t-y = 1 s ta = t-y. Right panel: same as the left panel except 
that ta = i^/lOO. Note that by decreasing the amount of time electrons have to radiate away their energy 
before being re-accelerated (ta) increases the max(no) roughly as t~^. The uq upper limit decreases when 
any of the GRB parameters {v^, f^, t^) is increased; no~f~^. 

understand. Frequent re-acceleration of charge particles makes it possible to have larger magnetic 
field while keeping z^c~100 keV. This decreases the number of particles required to produce the 
observed flux and that in turn makes it possible to have a smaller R^. 

We conclude that the synchrotron process in a shock heated medium cannot account for the 

p-i 

prompt 7-ray emission of long-duration GRBs with low energy spectrum (x. u 2 . However, 
synchrotron solutions appear to be viable when ta <^ t^, i.e. when electrons are accelerated repeat- 
edly, as might occur when magnetic field is dissipated and the energy is deposited in e^. 

3.2. Synchrotron solution when the low energy spectrum is us 

This is a special case of a = —{p— 1)/2 analyzed in the previous subsection (§3.1) when ~ jc] 
the solutions are a subset of those found in §3.1. The analytical solutions for this case, obtained 
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by substituting ,^ 



(see eq. [7|), fj. 



■^75! 



and f-y ~ fup, into equations 



are 



3.3 X l{)-^^v:;,^ f\t:^~Ha^Yl{l + Y)-^2{l + z 



3 

16 I 



p 



p 



1 

' 8 
''L2S 



" 16 



" 16 



Gauss 



^75 



1 i 

2d* 
J^28 



P 



P 



Substituting = ImJy, z/^^ 
these equations, we find 



1, tj = 0.1s (the observed parameters for GRB-\\), & ta 
r ~ 2.5 X 10^F"^(1 + 

7i ~ 2 X io'^yi^(i + y)3 

B 18Fi^(l + y)"t Gauss 

r ~ 6.3 X io"^°yi(i + 



(45) 
(46) 
(47) 
(48) 



ty, in 



(49) 
(50) 
(51) 
(52) 



and indeed, the solutions are a subset of the a = —{p — l)/2 solution space - these have smaller r 
and larger F & 7j. The distance of the source from the center of the explosion is: 



6.0 X 10^^ ul f^t'^tiY-l {I + Yf2 {I + z) 



P 



P 



cm 



(53) 



or Ry ~ 2 X 10 (1 + Y)2 cm for GRB-[\. This case has the same problems as a = —{p — l)/2 

case discussed in §3.1 i.e., large R-y and T, and requiring extremely small external density in order 
that Ry^Rd- Also, the conclusions drawn in §3.1 regarding ta/t^ <^ 1 offering a way out of this 
problem apply here as well. 

The numerical calculation of the hypersurface in 5-D parameter space allowed by GRB obser- 
vations - i^^, f-ykt^- for GRB-[\ finds 7^>10^ r>103, 200<7i/r<700, source radius (i?^) lO^^-lO^*^ 
cm, and B between 1 and lO^Gauss for the entire solution space (see fig. H]) - which is in very good 
agreement with analytical estimates. For a wide range of values for the three observable parameters 
we find the GRB source to be located between ~ 10^^ cm & 10^^ cm, 7j-Sj3 x 10^, and F>10^ (fig. 
U]). In order that R^/R^ < 1, the density of the surrounding medium (no) is required to be less 
than ~ 0.1 cm^^ which is much smaller than the value inferred from late time afterglow modeling 
for long duration GRBs. The density requirement is relaxed if ta <C (see fig. E]). 

In conclusion, the synchrotron process, in a shock heated medium, has serious problems ac- 
counting for prompt 7-ray emission for those bursts that have spectrum below the peak frequency 
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(u^) scaling as /jy oc i^^/^ or z^^^P^^)/^. A possible resolution is provided if electrons are more or 
less continuously accelerated while they are radiating 7-ray photons during the entire time pe- 
riod of a spike in the observed GRB lightcurve; in other words ta <C t^- It should be pointed 
that ta ~ ^7 in shocks whereas continuous acceleration might be possible in regions of magnetic 
reconnection / dissipation. 



3.3. Synchrotron solution when the low energy spectrum is u 

Substituting /j^p = (i/j^/i/cg)^ , u^^ = fj^, and ^ from equation ([7]) for the case where Uc < Vi, 
into eqs. [20H23] we find the allowed part of the 5-D parameter space when the spectrum below I'^y 
is oc 

T^w^4ff^tf'4Ay~^a+y)H^ + ^)'^di^^Af; (54) 

7 1111-, , -,1 1 

1 11 
^2p 



7, « 4.7 X 10^</7 'H^u-%tlYTE(l + y)4(l + z)4ci. M2 « (55) 
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Fig. 5. — Left panel: the maximum density of the circum-bust- medium so that < Rd for synchrotron 
sohitions with a = 1/3 and for a burst with i^-y — 100 keV, f-y = I mJy, t-y — I s &z, ta — ty. Right panel: is 
same as the left panel except that ta = i-y/100; note that by decreasing the amount of time electrons have 
to radiate away their energy before being re-accelerated (ta) increases the upper limit for uq roughly as t~^. 
The no upper limit is weakly dependent on v^, and decreases with increasing t^ & being most sensitive 
to - max(no)S/^. 



S«4.0i^^5«/^ ^H^Uc'Ha'^Yie (i + Y)—^{1 + z)-^d^lA^p^<^ Gauss 



r « 3.3 X 10 



-10, 



3 1 
8 ^8 



^75 



2p 



(56) 
(57) 



where 



2p 



p — 1 



(58) 



and taking < 10 (in 



p-2 

For GiJS-li ~ ty = 0.1s, Vy = lOOkeV, = ImJy t ta = 
agreement with the numerical calculation) we find r<720y~i6, 7j<6.7 x 10^, & B^590{1 + 
Gauss. In contrast to the previous two cases considered in §3.1 & §3.2, F < 1000 in this case. We 
find the 7-ray source distance, i?-y, to be 



Ry ^ 6.0 X 10 



16 



^« v^,tl y- 8 (1 + y) 2 (1 + z)-^dl^^Al^ cm 



(59) 
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which is cm for GRB-\\. We compare this radius to the deceleration 

radius, Rd, which is obtained from eq. (|26p and is given by 



2.6 X lO^*^ E^lgTiQ 3 1/^^8 H*fVc^^Ha''Ys{l + Y)-6{l + z)-(.d^^^A^^ cm s = 

(60) 

[1.8 X 1013 E^:,A-^v~'^ f~i^v~S~a~'Yl{l + y)-^(l + z)-2dllA~} cm s = 2 



and for the ratio Ry/Rd < 1, we find that no<8Er,3Y2(l + y)-2 cm'^ and A* < 0.014£;53y 3 (1 + 
y)~i if i^cg ~ 0.1; the limits on no & ^* are much higher for <C 0.1 and poses no problem for 
synchrotron solutions in a shock heated source. 

If the synchrotron solutions were to arise in a shock heated medium, we can calculate the LF 
of the shock front wrt the unshocked fluid, Tgh, using equation [TTl 

r,;, ^ 266-Vf /7^44iaV^(l + Y)Hl + z)UllAf^. (61) 

or r<j/i<32yT6 (1 + Y)i for GRB-\\ assuming that electrons receive half of the shock energy and 
that there are no pairs. Note that as long as v-^^ ~ 1, Tsh is pretty high (~ 20), and it is 
insensitive to z^cs (and the other quantities). In order to produce ~ 20 in internal shocks, we 



need the relative L F of the two colliding shells Tj-ei ~ 2r^^(ni/n2)2 ( Panaitescu &: Kumar 



2004 



Sari &: PiranI Il995l ). where ni and n2 are the comoving densities of the two colliding shells (see 
appendix A for a discussion of how we calculate Trei)- For Trei^5, so that the ratio of the LFs of 
the colliding shells is not larger than 10, ni/n2^ ~ 10^^ is required (see appendix A); T^ei > 5 is 
an unlikely situation to be realized in nature. 

We calculate the total energy in electrons {E±) and magnetic field (Eb) to determine the 
efficiency for synchrotron radiation — if there is a lot more energy in magnetic field than that in 
the electrons, the efficiency for 7-ray radiation would be small. The magnetic and electron energies 
for the case of z^"^/^ spectrum are obtained from equation ()10p and the solutions for B, T, 7j and 
R-y derived above, and are given by - 

Eb = B^R^Q « 6.0 x 10^° v^J^t^Y-^{l + z)-^dl^^A2peTgs, (62) 

and 

E± = N{p - l)meC^7iT/{p - 2) « 8.7 x 10^° ^^75/7*^(1 + y)(l + z) 4^8 ergs. (63) 
The ratio Eb/E± is 

fB/ke ^ 1^ « 0.68t^t-iy-i(l + Y)-\ (64) 

For ta ^ the above expression for the ratio fB/ke would need to be modified to include the 
total energy deposited in e^s as a result of multiple acceleration episodes during a GRB pulse time 
period of t-y; for t-y/ta ~ 1 the expression for fB/ke reduces to the familiar form that depends only 
on the Compton-y. For y ^ 1 most of the energy is in the magnetic field and for these solutions 
the radiative efficiency to produce a GRB is very small. 
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Fig. 6. — Synchrotron solution space when the spectrum below the peak of vfi, is fi, oc i/^^/^ for v < Ury. 
7,;, r, Y, and fs/ke S'^e plotted against the distance of the source from the center of explosion (R-y). The 
solution spaces for several different values of v-y, and GRB pulse duration (t-y) are shown in the four 
panels. Please see the Figure [2] caption for details. 

We numerically search the 5-D parameter space subject to the three observational constraints 
(v^, t^) and find solutions with 10^'^ cm < < 10^ cm, 10 Gauss < B < lO^Gauss, 100 < 
7i < 3 X 10^, 80 < r < 3000, and 10~^ <Y<7 (see fig. E]) — all in good agreement with analytical 
estimates presented above. We find 2 < Fgh < 100 and 10^^ < fs/ke < 10^- So it would seem that 
we have solutions with fB/ke ~ 1 ^sh of order a few - however it turns out that for /B/fce~10, 
Tsh^5. For Tsh^5, 10~^ < ^ < 0.1, and the ratio of the LFs of two colliding shells, Tj-ei, to produce 
this Tsh is greater than 20 (see fig. [7]) — fluctuations in the LF of the outflow with Trei on the 
order of a few are typically expected in internal shocks. 

Numerical solutions for the allowed part of the 5-D space for a range of observable parameters 
are shown in Figure [H An increase in z^^ leads to a slight increase of 7i and i?^ whereas F is quite 
insensitive to it. These behaviors are consistent with our analytical calculations (eqs. [53H57p . The 
decrease of fs/ke with (fig. E]) is due to an increase of Y . An increase of has little effect on 
7j (for the allowed solution space), F & y increase a little, and fs/ke decreases; these parameters 
have a very weak dependence on /-y (see eqs. [5lH57|) . And finally, when t.y is increased, 7^ and 
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20 keV 




Fig. 7. — This figure, unlike all the previous ones, shows model dependent results. For each point in the 
allowed region of the 5-D parameter space, corresponding to synchrotron solutions for a = —1/2 and for 
given fj & tj, we calculate parameters for the popular internal shock model for GRBs (see appendix 
A for details). The left panel shows the LFs Fi & F2 of the two colliding shells in the internal shock 
model as a function of i?-y, the distance from the center where the shells collide. The right panel shows 
Trei = rir2 (1 — V1V2) vs. fB/ke (the ratio of energy in magnetic field and electrons). Different v^, and 
cases are displayed as described in Figure [2l Note that F^e; is correlated with fB/ke', for high Treh fn/ke 
is small (< 1), and for low Trei, fB/ke ^ 1- Low Trei & fB/ke S> 1 solutions might correspond to a highly 
magnetized outflow. 



Y increase, and T and fB/ke decrease. This is again in agreement with the analytical estimates - 

- -- -A. 

7j oc . The large decrease in T with is due to an increase of Y with - T cct'^^Y le . 

We have looked at the variation of Tgh, ^rei and fs/ke with i/^, Sz t^. The results are shown 
in Figure [71 The minimum value of T^ei has a weak dependence on i^^, f^, and t^; Trei is between 
5 and 20 for fs/ke^^^'i ^rel ~ 5 solutions are only present when i^<0.01s, /-^^O.l mJy, or h'^ < 20 
keV (note that we only alter one of the three at a time, i.e. for ~ 0.01s, ~ 100 keV and /-^ ~ 1 
mJy). Trei may be small enough, then, that the synchrotron mechanism can produce GRBs with 
u~2 spectra if it has very short pulse duration, small peak frequency, or small flux. 
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3.3.1. X-ray flux during the GRB when a = 

So far, we have only considered prompt 7-ray emission due to the synchrotron process. We now 
calculate emission in other wavelengths, particularly the X-ray and optical, that should accompany 
7-ray photons. In this subsection, and in §3.3.2, we relate the solutions we found in §3.3 to the 
internal shock model for GRBs (Rees & Meszaros, 1994; see Piran, 1999, for complete references) 
according to which shells of material ejected in the explosion undergo collisions and the resulting 
shocks convert part of the kinetic energy of the outflow to radiation. Throughout this section 
we assume that the 7-ray emission is produced in shell '1' which is taken to be the faster of the 
two shells. The results are essentially identical if we assume that the GRB is produced in the 
outer, slower, shell, which we shall refer to as shell '2'. The x-ray flux from shell '1', however, is 
independent of the internal shock model, and is expected to accompany the prompt synchrotron 
V 2 7-ray emission. 

The x-ray and optical flux from shell '1' lie on the z^~^/^ extrapolation of the 7-ray flux and 
therefore it is straightforward to calculate these using and the information that i^c~2eV & fa<5eV 
for the entire solution sub-space of the 5-D parameter space. The calculation of emission from shell 
'2' is more involved and also a bit uncertain. We provide here (and in §3.3.2) a lower limit to the x- 
ray and optical flux from shell '2' for each point in the 5-D space that satisfies the three observational 
constraints (z^^, t^). The calculation of flux from shell '2' requires the knowledge of the LF of 
the shock front moving into this shell as well as the ratio of densities {n\l 712). The calculation for 
these quantities is described in appendix A. The synchrotron injection frequency in shell '2', i/jj, is 
smaller than that in shell 1 by a factor of (Fsj — 1)^ / (T^^ — 1)^; V & are shock front LFs into 
shell '1' & '2' wrt unshocked gas in shells '1' & '2' respectively. This factor is approximately ~ ^ 
for Fsj,Fs2 ^ 1, but the approximation breaks down for ^~10^^ (see Figure [T5]) since the shock in 
shell '2' becomes mildly relativistic (we note that this approximation is not used in our numerical 
calculations). The peak flux of the synchrotron spectrum at v = min ( z^jj, i^cj) in shell '2', fup^, is 

larger than f^^_^ by a factor of (^^^^ ■ The magnetic field is assumed to be the same in the two 

shells, and therefore the difference between the cooling frequencies in the shells is due to different 

^-parameters; since synchrotron dominates over SSC here by design, the difference between l'c2 

and ends up being very small. The shell '2' synchrotron self absorption frequency, t'aj is larger 

/ N 0.2 

than that in shell 1 by a factor of ( ^ j , or a factor of a few. 

The 1 keV synchrotron and SSC flux from shells '1' and '2' for GRB-\] is shown in the top left 
panel of Figure [8l The shell '2' synchrotron flux contributes the most to the prompt x-ray flux, 
and the shell 1 synchrotron flux contributes a slightly smaller amount. SSC flux from either shell 
is negligible. There is a weak dependence of x-ray flux on F^e/- The shell '1' synchrotron flux at 1 
keV is about 10 mJy. This is simply the extrapolation of the 100 keV flux back to 1 keV, with a 
spectral index of — ^ or (j^q) ^ ImJy ~ lOmJy for GRB-\\. 

Shell '2' synchrotron x-ray flux ranges from 1 to 100 mJy. We expect the flux from shell '2' at 
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Fig. 8. — Prompt x-ray [1 keV] and optical [2 eV] flux associated with points in the solution sub-space of 5-D 
parameter space when cast in terms of the internal shock model or two colliding shells; the solutions are for 
synchrotron radiation with low energy spectrum fi, (xv~^ . Top left: contributions of SSC and synchrotron 
to x-ray flux at 1 keV from shells '1' and '2' (see legend) for GRB-\\. Top right: sum of all contributions is 
shown for the cases that have been described in Figure [2l Bottom left: contributions of synchrotron emission 
from shells '1' and '2' to the optical flux; IC contributions are negligible from either shell. Bottom right: 
sum of all contributions in the optical, for the same cases as the top right panel. The results shown here 
assume that the 7-ray flux is produced in shell '1'; production of the GRB in shell '2' does not change the 
results significantly. 



/lkeV\ 2 / zyj2 \-5 „ /ni\4 /ikeV\ 4 



1 keV to be 

using the above arguments and assuming that < 1 keV (valid when ^<0.01 - satisfied by 
roughly half of the solution points in the 5-D space), ~ t'cj, and p = 2.5. Prom this equation, 
we estimate that the flux at 1 keV should be between ~ 20 and 300 mJy (for fi^ = 1 mJy and 
~ 100 keV) - in agreement with the numerical results in Figure El 

The shell '2' SSC flux at 1 keV ranges from 10^^ to almost 10 mJy. The SSC peak frequency 
(~ 1^42 ranges from about 100 eV to very high values, so over a large part of the solution space. 
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the expected SSC flux at 1 keV in comparison to the synchrotron 100 keV flux from shell '1' is 
(assuming that > i'c2 ) 

A fvn V7c2^c2/ VlOOkeV/ V"i/ 

which, after substituting in the solutions for r and 7c, is 

, f rio\ 1- 11 _i _i _i 15 1 1 1 _i5 

fl « 1.5 X 10-^ (^-^j v^-^ f-HY'yc,Ha~'Y-.{l + Yr-^{l + z)"dl^^A^;'^ mJy; (67) 

with y ~ 1 and 1 < '/12/ni < 10^, the range for the x-ray flux obtained from the above equation is 
in agreement with the numerical solutions. 

The sum of synchrotron and SSC contributions to flux at 1 keV from both shells for various 
values of z^^, and are shown in the upper right hand panel of Figure [HI The 1 keV flux ranges 

from 0.1 to a few thousand mJy, and is most sensitive to and - in agreement with equation 

5 

([65]) . oc /^i^-f , since synchrotron emission dominates. 

The early 0.2-10 keV x-ray flux as observed by the Swift x-ray telescope ranges from 10"^'^ 
to 10~^ ergs cm~^ s~^l£| which corresponds to 1 keV flux of about 10~^ to a few mJy (assuming 

in the x-ray band). These observations are made at roughly 100s after the GRB trigger. The 
x-ray light curve from the 7-ray source should peak at about the same time as the GRB light curve. 
After the peak, assuming that the outflow opening angle is greater than l/F, t he emission should be 



dominated by off-axis emission and the light curve should fall off as t ^"''^ ( Kumar &: Panaitescu 



2000), which in this case is since fy oc v'^l'^ . Extrapolating the observed 1 keV flux of 

~ 10~'^-1 mJy backwards in time from 100s to 10s, we find the x-ray flux during the GRB to be 
consistent with values shown in the upper right panel of Figure [H 



3.3.2. Prompt optical emission when the GRB index a = 

In the bottom left panel of Figure [HI the R band (2 eV) flux from shell '1' and '2' are plotted 
against Trei- Optical flux from the 7-ray source can be pretty bright during the burst for these 
solutions ranging from 10"^ and 100 mJy (24th to 11th magnitude in the R band) . The synchrotron 
flux is smaller for smaller Trei solutions. The SSC makes negligible contribution to the optical flux 
compared to the synchrotron process, because fa 7c is well above the optical. 

If we extrapolate the shell '1' 100 keV flux back to 2 eV using the spectral index q = — ^, we 

1/2 

expect fiR ~ 2251^75 ^ 225mJy for Gi^iJ-t] whereas for most of the solution space the optical flux 
for shell '1' falls below 10 mJy (fig. [8] bottom left panel) — this is because the synchrotron self 



^urlhttp: / /swift .gsfc.nasa.gov / cgi-bin / swift / grb_table /grb_table.py 
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absorption frequency is larger than the R-band frequency by a factor ~ 10 or more. The range of 

shell '2' R-band flux is higher than shell '1' by a factor of ~ ~ 30. 

ji^Pi y '''1 

In the bottom right panel of Figure El we show the affect of varying i^^, and t-y on the 
prompt optical flux; the total flux - obtained by adding the contributions for the two colliding 
shells - ranges from 10~^ to 10^ mJy, or R magnitude of 26 to 9. The optical flux increases with 
u^, f^, and ~ longer GRB pulses with higher peak frequency and/or flux tend to be brighter in 
the optical band. Synchrotron self absorption is larger at smaller Trei, and that makes the optical 
flux smallest at the minimum of Trei- There are Tr(,i ~ 5 solutions with fs/ke ~ 1) that have small 
enough optical flux to be in accord with the observed upper limits, especially for smaller v^, f^, 
and t-f. Although fa > 2 eV for much of the solution space, the optical light curve should peak at 
the same time as the GRB light curve, since f c < 2 eV. After the peak, the light curve should fall 
off - if dominated by off-axis emission - as t~'^'^. 



3.4. Synchrotron solutions for fyccv 2 

If the observed spectral index is steep and consistent with v 2 and no break is detected in the 

_ 1 (p-i) 

observed energy band of 15-150 keV, for instance, then this is a special case of either 2 or 2 
low energy spectrum discussed in §3.1 & 3.3 - with 1/^5 = max(z^^5, t'c5)^0.15. The allowed solution 
space for this situation should be close to the = 20 keV case in Figures [2] and Specifically, 
7i<103, F > 100, and 2<Trei'^2m. 



4. Synchrotron-self-Compton — SSC — solutions 

In this section, we present solutions for the prompt 7-ray emission to be produced via the 
synchrotron-self-inverse-Compton radiation or the SSC process. The basic approach is same as 
in section 3. We determine the hypersurface in the 5-D parameter space {jijT, B, N,t) that has 
SSC emission consistent with the three observational constraints ly^, t^. Since different cases 
of low energy spectral index have different ordering for the characteristic synchrotron frequencies 
(vijUcUa), we do not consider a general SSC solution, but describe analytical and numerical solu- 
tions for the positive low energy spectral index case i.e., fi, oc z^" with q > for 1/ < i'^, and the 
negative index case i.e., a < 0, separately in several subsections below. 



4.1. SSC solutions: Positive low energy spectral index 

This section is broken up in two subsections. One dealing with the special case of f,y oc z/^/^ is 
discussed below. All the other cases of a > 1/3 are discussed in §4.1.2. 
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4.1.1. SSC solutions: a ~ +1/3 

The SSC spectrum {vfl^) peaks at ~ 4max(fi, fc) max(7j, 7c)^; where Vi & Uc are the 
injection and coohng frequencies of the underlying synchrotron radiation, 7j is the minimum LF 
of electrons in the source comoving frame and 7c is the LF of electrons that cool on time scale ta 
available since last accelerated. For the spectrum below to be ~ v'^^'^, we must have 7^ ~ 7c, 
and in that case ~ ^I'i'yf . The IC flux at is /-y ~ /^pT {f^^^ is the synchrotron flux at Vi). The 
equations for pulse duration and Compton-Y are same as in §3. The equations for z/j ~ Vc, the 
peak IC frequency z/^, the IC flux at z/^, and the Compton Y-parameter are given below - 

B'^T-ii w 7.7 X 10^(1 + z)t~^{l + Y)-^ (68) 

B-ffV « 2.3 X lQ^^u^.^{l + z) (69) 

BFV^ « 1.6 X 10^7*^^(1 + z)dl^^ (70) 
2 /^P- 1 



(71) 



We first eliminate r from equation ([70]) using ([7T]) to get 

2 

Next, we divide equations ([55]) and to eliminate F: 



i?F57r' « 2.8 X 106/7*^^(1 + z)d%J-^ ( ^ ) . (72) 



S7-' ^ 3.3 X 10-^t-i(l + y)-V-\ (73) 



75 

divide equations ([69]) and ([72]) to eliminate i?: 

-2 



p — 1 
p-2 



^fr-'^8.2xlO\j-H'^dllY'['-^] , (74) 



and combine equations ([68]) and ([72]) to obtain: 



Equations ([73]) and ([75]) give 

2 

F73 « 1.4 X 10«i.^,i|(l + Y)lfiH-hl + z)ldlY~l ' , (76) 

and substituting this into equation ([74p . we find the solution for 7^ to be 

7i « 844/7"^*?4l^^(l + Y)U-^^^1 + z)U2„'^- (77) 
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Note that the electron LF 7, has a very weak dependence on the observed quantities as well as 
the Compton-y parameter, and therefore 7^ ~ 80 for the entire SSC solution space. By plugging 
equation ([77|) back into equations ([73]) , ([71|) , and ([7T]) , we find the remaining parameters 

B w 2001/^3^7 ^4*a^^^(l + Yy^'dllil + z)^A~^ Gauss (78) 
r « 24o4 # t-^HiY-T-sil + Y)! d (1 + ^) i A j (79) 

r « 1.1 X w-\^/^ftH-h-hl{i + y)-i4s(i + ^)'^^27- (80) 

All of these parameters are weakly dependent on the three observable quantities viz. u^^, &; 
try. Substituting the observed values for GRB-\\, i.e. z^-y, = 1, = 1 mJy, and ~ = 0.1 
s, z = 1, k taking p=3.2, we find that 7^ ~ 58^1^(1 + F)^, S ~ 1.1 x lO^Y^l + y)~3Gauss, 
r ~ 110y~^(l + y)i, and r ~ 1.3 X lO^^yt (1 + F)"! . 

The distance of the 7-ray source from the center of the explosion, R^, is given by 

R^^3.3x 10^^vlf^4tiY-^l + Y)l{l + z)-Ul^^Aj^cm, (81) 

or Rj ~ 1.3 X 10 "'y^g (1 + y)9cm for GRB-\\. This distance is smaller than the deceleration radius 
for a homogeneous or a wind external medium, unlike the situation when 7-rays are produced via 
the synchrotron process (see §3). 

One constraint that we have not yet considered is that the SSC self absorption frequency, 
ulf ~ 4fa7?, must be smaller than ~20 keV otherwise the low energy spectral index, obtained 
by Band function fit to the BATSE or Swift /BAT data, would be steeper than a = 1/3 we are 
considering in this subsection. The expression for v^^, valid for Va < '^i,t^c, is 

1.7 X io-^s?rt / /;-r' , .... 
irT7) — 1^2^ I 

where fp = ^/^q^ Bt / aTm^c? is the comoving synchrotron flux. Substituting for 7^, T, i? &: r using 
equations iTTtiSOl we find 

1^:^-2.2 X lO^i^f /|t7t;^yi^(l + y)-ii(l + z)-i4«^2,^eV (83) 

which is very insensitive to all of the observed quantities and for a wide range of observables 
f*'^ ~ 100 keV which is too large to produce an SSC spectrum with fif oc z/s below the peak. 

We now try relaxing one of the constraints we had imposed to simplify the analytical calculation 
i.e., 7j ~ 7c. This approximation was guided by the observational result that the observed spectrum 
for > z^^ is almost always ~ for GRBs with a ~ 1/3. This result suggests 7^ ~ 7c provided 
that p ~ 3. However, if the electron distribution is steeper, p ~ 4, then Uc can be much greater than 
z/j, and a high energy spectrum of would be consistent with observations. We now 
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investigate this possibility and determine if letting > t^i would allow for a smaller i'^ and hence 
a = 1/3 solutions. Note that the opposite arrangement of frequencies (fc < fi) is uninteresting, 
since the low energy index is -1/2 in this case. 

For i/c > equation ()68p is modified to read 

where rji = jc/li- We also need to use the appropriate expression for Y when 7c » ji and p > 3 
(see eq. [6]) 

4 jp-l) 2 /o.x 

^=3(,-2)(,-3)-^- 

which apart from a factor (p — 3) is same as equation ()7ip . We solve the above two equations 
_j_ 

together with eqs. ([70|) & ([7T]) to find that z^*^ oc r]- ; so i'*'^ does not decrease by much even if we 
take 7c to be larger than 7^ by many orders of magnitude, and therefore there are no SSC solutions 
with low energy spectrum of v^^^ between ~ 15 keV and 200 keV. 

The above analytical calculation is based on a number of approximations for the SSC spectrum 
and flux. We check the validity of analytical results using numerical calculations and by searching 
the 5-D parameter space for SSC solutions with low energy spectral index a ~ 1/3. It turns out 
that numerically also we find no solutions — is indeed too high to produce a GRB with a low 
energy spectrum of i^^^^ in the ~15-200 keV band. 

The only other possibility is that the a ~ 1/3 index is transitory, i.e. the spectrum is changing 
continuously from a ~ 1 at ~ 20 keV to a ~ — 1 at z^ = i^-y and that a ~ 1/3 is realized at some 
intermediate frequency. This might, however, pose a problem for those GRBs with z/^>100keV, 
since the spectrum would be steeper than v^^^ near 15 keV and therefore a Band function fit to 
the spectrum will yield a > 1/3. 



4-1-2. SSC solutions: Spectral index 1/3 < a < 1 

The analytical solution for this case is similar to the SSC q = 1/3 case analyzed in §4.1.1. We 
take 7i ~ 7c in order that the high energy spectrum is oc v~^l'^ ~ v^^'^ . The equations we solve 
are for i/^, Compton-Y, and ^c- 

B^T-ii « 7.7 X 10^(1 + z)t-^{\ + Y)-^ (86) 
B^tV ^ 2.3 X 10i2z.^^(l + z)^-2 (87) 

i?rV2 « 1.6 X ioV7t7'(i + ^)4,s^a 



2 / P - 1 



-1 
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where rja = Ja/li^'^ since Ua^i'c, i^i- The above equations are solved in the exact same way that we 
solved them in section 4.1.1, and we find that the solutions are: 



7i « 84</^ ^«^n„«yTB (1 + Y)-9d^^\^ (1 + z)-9A^;'rja 



-(p+18) 



36 



1 6-p 



B^ 200i^jr/fy ^H^ta''Y6{l + Y)-^d^l{l + z)3A^p''ria^^ Gauss 



, 7 7p-18 
^ —36— 



8 P+18 
18 



(90) 
(91) 
(92) 
(93) 



r ^ 2A0iy4J^^H^ 'HiY--s{l + F)9(if^^(l + z)9A^;rja 

r « 1.1 X W-%,^ftH~Hahl{l + Y)-ldlJl + z)-^9A-;riP 

For p = 3.2, ta ~ t-y, ^a^^l, and parameters corresponding to GRB-\\ we find from the above 
equations that r>110y^i8 , 7j<100, and t>1.3 x 10^^y9(l + Y)~9 . These analytical results are 
roughly consistent with numerical determination of the allowed region in 5-D parameter space (see 
fig. [9]); we also find r/^ ~ 1 numerically. 

The distance of the 7-ray source from the center of the explosion is shown in fig. ([9|) for various 
GRB parameters and is greater than ~ lO^^'cm for y<10. The ratio of magnetic to electron energy 
is small - fs/ke < 0-1 for the entire solution space (fig. [9]) since Y>1. 
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Fig. 9. — The SSC solution space when the low energy spectrum is /^'^ oc i^", for < v^, with 1/3 < a < 1. 
The different cases in legend are as described in Figure [2l 



The SSC solutions we have found can be related to the internal shock model. The relative 
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LF of collision between shells - obtained from 7^ (see eq. [TT] and appendix A) - is found to be 
between 2 and 10, which is significantly less than what we were finding for synchrotron solutions. 
The LFs of shells before collision (assuming that 7-rays are produced in the inner, faster, shell) 
is 300<ri<5000 and 100<r2^1000; the ratio ri/r2 > 2 and the efficiency for producing 7-rays is 
>10% for the allowed 5-D parameter space for SSC. The bulk LF of post-shock gas 100 < F < 1000 
is compatible with late time afterglow modeling. 

We now calculate the x-ray and optical emissions accompanying the 7-ray pulse. 

4.1.2a. X-ray emission for 1/3 < a < 1 SSC solutions 

The 1 keV prompt emission from SSC and synchrotron processes is shown in the top two panels 
of figure [TOl The contributions of SSC & synchrotron to 1 keV flux is shown separately in the top 
left panel for GRB-\\, and the sum of the two for a variety of GRB parameters can be found in the 
top right panel. 

The x-ray flux can be estimated analytically using the expression for synchrotron fiux fx = 
/t,p(lkeV/fj)~^/^, since i/j ~ i/c ~ Va\ fx can also be expressed as fx ~ /7T~^7~^(25)P/^f^g^, or in 
terms of observable parameters 

q X in3 p+2 34+p 2-p 2-p ig, 2_ £-2 2-p (p-2)(p+18) 

/. ~ /7 V' ta' Y- — {1 + y) — <3 (1 + z) — A^l^ i^a mJy, (94) 

16 2 

or fx ~ 9y~i5 (l+y)~ 15 mJy for GRB-\\ with 7]a > 1, p = 3.2 and ta ~ t^; this is roughly consistent 
with the numerically calculated fiux shown in fig. [101 

In the top right panel of Figure [T0| the total 1 keV prompt flux is shown for a number of 
different values of i^^, f-y, and t^. The dependence of the x-ray flux on these quantities agree with 

13 31 _J_ 

equation (j94|) . which gives fx'^'^-y^ f'y^' t-y for p = 3.2: an increase in 1^^ or leads to an increase 
of fx, and an increase in has little effect on the x-ray flux (fig. [ro]) . The 1 keV fiux for all of 
these cases ranges from 0.01 to 10^ mJy during the burst; the fiux at 100s, the time when the x-ray 
telescope aboard the Swift satellite starts looking at the burst, would be smaller by a factor of 
~ 10-10^ depending on GRB pulse duration. So, the x-ray fiux accompanying the 7-ray radiation, 
for the SSC model of GRBs, is consistent with the observed data. 

4.1.2b. Optical emission for 1/3 < a < 1 SSC solutions 

The prompt optical fiux accompanying 7-rays, in the SSC model, is shown in the bottom left 
panel of Figure [TOl Analytically we find the prompt R-band flux due to the synchrotron component 
underlying the SSC model to be: fji ~ f^p ~ f-yT'^, or 

1 17 1 2 g 2 _i 2 ^ 17p-18 

~ 9 X 1034/«t|i|y-9 (1 + Y)-9d^l{l + z)-9Alr,a mJy (95) 

which for GRB-\\ reduces to f^ ~(8xl03 mJy) Y-9{l + Y)9. There are, however, many numerical 
solutions corresponding to Y>10 for which //j ~ 70 mJy. It turns out that > 2 eV for the low 
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Fig. 10. — X-ray (1 keV) and optical (2eV) synchrotron and SSC flux accompanying the prompt GRB 
emission for SSC solutions with 1/3 < a < 1. Right panels: The sum of SSC and synchrotron contributions 
for several different sets of v^, & as described in Figure [8] Left two panels: the areas shaded in gray 
are the x-ray and optical flux for GRB-\\ (a GRB that has ~ lOOkeV, — ImJy Sz t-y = 0.1s) and with 

ta ~ t-y/lOO, and the unshaded areas are for GRB-\\ with ta ~ t^. 



radius solutions, by up to a factor of 3. If i/j < 2 eV< f^, then the expression for optical flux is 

_3 29 J7_ 7 11 y _T_ 7 11 (7p+126) 

fR^f>Qv^.^f^H),HlY--B{l+Y)-^d^ll{l + z)i>A^^lrja mJy (96) 

giving fR<7Y~i^{l + Y)l mJy for GRB-[\. 

The optical flux, obtained by numerical calculations, is shown in the bottom right panel of 
Figure [TOl for several sets of (z^^, f-y, t-y). The results are consistent with the dependences found 
in equations (f95]) or ([96]) when > 2eV. The reason that the self absorbed /r increases with z^^ 
numerically while equation (|96p shows a decrease is that rja ^ v-y (confirmed numerically), and 
the huge dependence of /r oc t/^'^ gives positive dependence of /r ~ u}^'^. The dependence of optical 
flux on the duration of the GRB pulse is due to the fact that longer pulses have larger R.y and 
Ua < 2eV. The range of optical flux for SSC solutions is between 0.01 mJy to a few times 10'^ mJy 
(i? magnitude from 21- to 6-mag). There is an approximately linear relationship between Jr and 
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R^. Solutions with i?-magnitude of above 9-mag (1 Jy) are most likely ruled out. In particular, this 
rules out the 1 s pulse duration solutions that have ii^>2 x 10^^ cm. If the pulse width were 10 s, the 
SSC solutions would have prompt optical of between 1 and 4 Jy, or i? ~ 7 mag, which is too bright 
to have been missed ir i optical follow up o bservations. We note that if GRB dissipation radius is 
~ lO^^cm as found in iKumar et al.l ( 20071 ). then bright optical flux of i? > 9 mag is expected in 
every GRB produced via SSC. Since this bright optical emission is not seen, this may pose major 
problems for the SSC process to produce GRBs with positive a. 



If, however, electrons are accelerated multiple times during the course of a 7-ray pulse in 
the GRB light-curve, i.e. ta <C t^,, the optical flux can be reduced significantly. The dependence 

4(p+3) 

[. [95]) itself does not reduce fn by much, but since Va oc ta ^'^^ ~ ta^'^ as well, a 
smaller ta gives a larger Va and that reduces the optical flux by an additional factor of ~ t~^'^ 
and results in fR^tU'^; for ta = t-^/lOO, the optical flux is reduced by a factor of ~ 10^ compared 
with the case where ta ^ t^, in agreement with the numerical results found in the lower left hand 
panel of Figure [TOj We've numerically searched the whole range of observable parameters and find 
that the scaling fu^tj/^ is valid through the entire solution space, and even the highest optical flux 
levels of 10 Jy (for the t^ ~ Is case) is reduced below 0.1 Jy, or i? > 11 mag, if ta<^/10^. Multiple 
acceleration episodes for electrons is, therefore, a possible way of reducing the excessive optical flux 
that otherwise necessarily accompanies SSC solutions for 7-ray emission. 

The optical flux should peak at the same time the GRB light curve peaks, since Vc is on 
the order of 2 eV. The temporal decay of optical flux in this case is dictated by the curvature 
or the off-axis emission (Kumar & Panaitescu, 2000) — the optical light curve should fall off as 
i-2-p/2 ^ ^-3-6 iQfr^g jg below the R band. At first glance, it might seem that if Ua > 2 eV, 
and the synchrotron spectrum oc i'^^'^ or z/^ in the optical band, the optical LC would be flat ~ t^ 
or even rise as t^^'^. This behavior, however, lasts for a very short time since (fa/2eV);S a few, and 
Ua oc t~^ for off-axis emission; once Va drops below 2eV the optical flux would start falling off as 
^-3.6 rpj^g upper limit of y ~ 18.5-mag (0.2 mJy) at 100s for many Swift detected bursts (e.g. 



Roming et al.ll2006l ) is a lot smaller than the flux expected during the burst for the SSC model. 
If the pulse occurs at ~ Is post-trigger then the optical flux at 100s would be smaller than the 
prompt optical flux by a factor ~ 10^ and that is quite consistent with observational upper limits 
for almost the entire solution space for the SSC-model. 

We emphasize that a bright optical flash (i2<14-mag) concurrent with 7-ray emission is a 



optical radiation has been reported for a few bursts wit 


1 positive a - e.g. 990123 ( 


Akerlof et al. 


1999; 


Briggs et al. 


1999). 061007 ( 


Golenetskii. S., et al. 


2006; 


Yamaoka. K.. et al. 


2006 


), and the 


second emission episodes of 050401 


Golenetskii. S.. et al. 


2005) and 050820A ( 


Cummings. J., et al. 



2OO5I ) - however, if future observations fail to detect prompt optical with i?<14-mag then that will 
suggest that one of the assumptions of the model developed in this work has to be abandoned - 
the most likely possibility, in our view, is to discard the assumption that ta ~ t^ and that suggests 
that 7-rays are not generated in a shock heated medium. 
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4.2. SSC solutions for negative low energy spectral index 

In this section we consider synchrotron-self-Compton solutions when the spectrum below i/^, 
the peak of lyfl^, is fl^ oc with a < 0. There are two different class of solutions in this case 
- those with the underlying seed synchrotron spectrum having Uc < fi and vice versa. We treat 
the two cases separately analytically, but plot the numerical solutions for both cases together in 
Figure [TTJ We use one vital piece of information gained from the numerical solutions to simplify 
our analytical calculations: the synchrotron self absorption frequency, fa, is larger than max(z/j, v^), 
and therefore ~ 4fa max(7j, 7^)^; note that even though > max(fj,i/c), the spectral index 
below is negative down to the frequency ~ min(7j, 7c)^ < 10 keV. 



4-2.1. f c < Vi case 

The equations that are solved for this case can be cast in a form very similar to the SSC 
a PS 1/3 case (considered in §4.1.1) by introducing two variables: r/j = 7i/7c and r]a = 'Ja/ji- The 
equations for Uc, v^, f-y, and Compton-y expressed in terms of r]i & r]a are: 

B^T^i ~ 7.7 X 1087?i(l + z)t-\l + y)-i (97) 

Bjfr ~ 2.3 X IQi^zv^, (1 + z)r,~^ (98) 

i?rV2 ^ 1.6 X W'fyt~\l + z)dl^y,rjf (99) 

T^f ~ 0.75Y Vi- (100) 

These equations are solved the same way as outlined in the previous section, and we find the 
solutions to be: 

7. ~ 84u}^i'f-'/^%/'X^'Y'^'Hl + Y)'/'dll/''{l + zfl^A-'J'^Prif^ (101) 
B ~ 2mu-^l'f-'''hyh-JI^Y'/\l + Y)-^l'dlll\l + z)"'A^"riri? Gauss (102) 
r ~ 2mu]l^ fy^H-y^h^fY-y^\l + Yf/^d'l^^{l + zf^A^^'r^Pri^ (103) 
r ~ 1.1 X f}^l'%yVl^Y^/\l + Y)-'''dZ{l + z)-''' A-^'\f r^a^ . (104) 



The dependence of B, T, ji h t on the observables is same as in equations iTTllSOl ~ the difference is 
in the dependence on r]a and r/i. The distance of the 7-ray source from the center of the explosion 

_ 4 7p-18 

is R-y oc t^r^ oc 7] ■ rja . We use the constraint that ~ 4fa7^ ~ lO^i^^giy"^ < 10 keV to infer 
that rii^3.2 for these solutions; we numerically find that T^j^^lO. 

Numerical calculations of the allowed region in the 5-D parameter space for GRB-\\ give 
7j<50, -B-^200 Gauss, r<300, and r>3 x 10^^. There is good agreement between the numerical and 
analytical 7^ and T solutions (see fig. [TT]) . although analytical and numerical solutions for B, t &i 
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can differ by a factor of a few due to the difference of a factor of a few in the analytical and 
numerical value of r]i. 



20 keV 

— — — 400 keV 
0.1 mjy 




[cm] [cm] 



Fig. 11. — Synchrotron-Self-Compton solutions when the spectral index, a, below (the peak of vfl'^) 
is less than 0. There are two branches of the solutions both of which are shown in the figure. For one of 
these branches Vc < Vi, and for the other Uc > i^i] 7-ray sources corresponding to the Vc < Vi branch lie at 
a smaller distance from the center of explosion {R^) than the other branch. Allowed regions for the 5-D 
parameter space is shown for a number of different sets of GRB observable parameters (see figure [2] caption 
for details) . The top right panel shows the LF of the 7-ray source (one of the five basic parameters we use 
to describe the source) as well as the Trei - the relative LF of collision of two shells obtained by mapping 
the 5-D parameter solution to internal shocks (see appendix A). 

In figure [TTl we show the numerical solutions for the allowed region in the 5-D space for the 
SSC-model for a wider range of observables. Note that results for both i^i > Uc I'i < Vc (to be 
considered analytically in the next subsection) are plotted together in fig. [TTl The allowed range 
for electron LF 7j is 3-200 (fig. [TT]) which is characteristic of mildly relativistic shocks. If we cast 
the 5-D parameter solutions in terms of colliding shells as described in appendix A, we find the 
relative LF of collision between shells to be less than a few. For F ~ 100 &: Trei ~ a few, there is 
little chance of an external forward-shock origin for these SSC photons. However, the 5-D solutions 
we find appear to be consistent with an internal shock; is smaller than the deceleration radius 
and fs/ke^^ for ttis entire SSC solution space. 
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The dependence of the solution space on z^^, and is in agreement with analytical estimates; 

1 __L 1 

for example, we have verified the analytical dependence, 7, oc 2^7/7 '^^t^ (for ta = ty), given in 

equation (jlOip . and T oc f^^t^ ^ (eq. I103p . We find that SSC can produce a wide variety of 
GRBs with the low energy spectral index a ~ —1/2 without requiring extreme parameter values. 

In order for these solutions to remain viable, the prompt optical and x-ray flux needs to be 
in accord with observations and/or upper limits. We next look at the prompt x-ray and optical 
emission from these solutions. 

4.2.1a. Prompt x-ray emission for SSC solutions with a < —1/2 

In the top two panels of Figure [I2l we've plotted the synchrotron and SSC contributions to 
the x-ray (1 keV) flux accompanying the 7-ray emission during the burst. In the top left panel, we 
see that the 1 keV flux ranges from 0.1 to over 100 mJy for GRB-\\ and much of it is due to the 
underlying synchrotron emission. 

There are at least a few solutions for each value of 1^7, f-y &i t.y we have considered with x-ray 
flux less than 10 mJy (fig. [T2|) with the exception of the f.y ~ 10 mJy case. The high end of this 
range of x-ray flux is above the value typically observed by Swift at 100 seconds following the burst 
(10~^ to a few mJy); but we know from early x-ray observations that this flux is initially falling off 
very steeply, ~ and therefore x-ray flux of ~ 10^ mJy during the burst, t7<10s, would be less 
than 1 mJy at 100s, or within the observed flux range of the x-ray telescope aboard Swift. 

We expect the x-ray flux to be 



1 



AkeV\-V.A-^ (105) 



(since synchrotron dominates). In terms of the observable quantities the flux is 

n X in3 p+2 34,+p 2-p 2-p ig , 2-p 2-p (p-2)(p+18) p-2 

f.^—^vy,' /7^« ty^'ta' (1 + y) — (i4(l + z) — A^p" T,^' mJy, (106) 

so for p = 2.5 and ta ~ ^7, fx oc f^'^t^ . This analytical formula is in good agreement with 
numerical solutions shown in fig. [T2j 

The upper limit to the synchrotron flux can be understood from the limit we place on the 
synchrotron 10 keV flux. We restrict synchrotron flux at 10 keV to be less than the 10 keV SSC 
flux, in order that the GRB spectrum is SSC dominated. The 1 keV synchrotron flux is then 
< /7(10keV/z^7)(P^-'^)/^ ~ 55 mJy for GRB\ in agreement with Figure [T2l 

4.2.1b. Prompt optical emission for SSC solutions with a < —1/2 



In the bottom two panels of Figure \T2\ we plot the optical emission accompanying the 7-ray 
radiation for the SSC solutions. In the optical (2 eV), we look only at the synchrotron emission. 
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Fig. 12. — X-ray (1 keV) and optical flux simultaneous with the GRB emission for the SSC solutions with 
a < are plotted against the 7-ray emission radius. The top left panel shows the synchrotron & SSC 
contributions to the 1 keV flux for GRB-\\, and the bottom left panel shows the optical (2 eV) flux associated 
with GiJS-tj-SSC solutions. The right panels show the 1 kcV flux (top right) and 2 eV flux (bottom right) 
by summing up contributions of synchrotron and SSC for various sets of observed parameters for GRBs as 
described in the caption for flg. [21 



since the SSC flux is highly suppressed as ^ 2eV. The optical flux is between 0.2 and 100 
mJy, or between ~ 18 and 11 mags for GRB-\\ (fig. [T2|). For the range of i^^, and ^ we 
have considered, we find the optical flux to be between 0.02 mJy and 50 Jy, or i? ~ 20 to 4 mags 
(bottom right panel of fig. [T2l) . Observational limits on R ban d flux, 100s after th e burst, are 
>18.5 mag, or ^0.2 mJy, for approximately half of Swift bursts (jRoming et al.1 l2006l ) . For small 
values of u^, f-y, and there is no problem satisfying this upper limit, but large v^, f-y, or might 
exceed the observed optical flux limit. Since the prompt optical flux falls off very rapidly with 
time as for t > tj, SSC solutions with hundreds of mJy optical flux during the burst are 

consistent with the upper limit of ~ IBmag at 100s. The solutions with optical flux greater than 
100 mJy or so, however, can be ruled out, since prompt optical flux this bright is very rare. 



Analytically, the R-band flux, dominated by synchrotron photons, is 



(107) 
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provided that Va > 2eV, and i^c < t'j < fa- In terms of the observed quantities the flux is given by 

_ 3 29 _7_ 7 j^j^ 7 _X 7 ii ^(p+lS) _7 

/r ~ 2 X 10=^ l/^5V7''*7'*al^"^(l + ^)^^Li' (1 + ^)^^2p^a ^^ ' ^JY- (108) 

_ 3 29 7 

If we assume that ~ we find oc i^-y ^ f^'^t^. This agrees with what we find numerically for 
and t^, but not v^. Numerically, increases with v^, while this expression shows a decrease. 
This difference is caused by the sensitivity of //j ~ r/~^ - numerically we find that rja oc 

5 

changing the above dependence on z^^ to be //j ~ z^^ , in accord with results shown in Figure [T2l 

11 7 

For GRB-\\, //j<120y^ is (1 + y)9 mJy. This estimation for fji is larger by a factor 10 than 
what we find numerically (see fig. I12p . This factor of 10 difference is due to the fact that 7c ~ 1 
for many of these solutions. When 7c < 2, we have a population of electrons that don't radiate 

synchrotron emission and we need to reduce the number of radiating particles. This is done by 

1 

using fc (7c — 1)^. Since fn oc fc , the correct value of //j, is a factor of 7c/(7c — 1) smaller than 
the crude analytical estimate above. The smallest value of 7c that we find numerically is 1.08, and 
therefore the analytical expression for overestimates the true flux, calculated numerically, by a 
factor of about ~ 14. This indeed reconciles the analytical and numerical results. 

The R-band flux increases with increasing t^. This is because increasing the pulse duration 
increases the radius at which the GRB emission is produced, and the synchrotron self absorption 
frequency is smaller at larger radii, and therefore brighter optical flux is observed. Thus, a prediction 
of the SSC model is that brighter optical flux accompanies wider GRB pulses — very spiky light 
curves (with short variability time scale) will have small optical flux that can escape detection, but 
those with wide pulses should have bright early optical afterglows. If a pulse duration were to be 
10 s, we should expect prompt optical flux of 100 mJy or larger (R-band magnitude smaller than 
ll*''). If this optical emission is not detected, it will point to one of the assumptions in our model 
for the SSC-emission being violated - most likely ta 9^ t^, i.e. electrons are not accelerated just 
once, but multiple times, during the course of a pulse duration of t-y. 



4-2.2. Vi < Uc case 

Analytically this case is very similar to the case of Vc < '^i discussed in §4.2.1. We require this 
time that 7j/7c = rji < 1 and r]a is still > 1. Since the ratio 7a/7c = Va^ii we will see that some of 
the solutions can be expressed in terms of r^a^i instead of using r]i and r]a separately. 

The equations that we are solving for, in this case are 

B^r^i ~ 7.7 X 10\i(l + z)t-^{l + y)-i (109) 

Bjfr ~ 2.3 X 10'^u^,{^ + z)^-\f (110) 

BT\^ ~ 1.6 X 10V7i7'(l + (111) 
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These equations are almost identical to the SSC z^a case of §4.1.1, with exception of the dependence 
on the variables rji and rja and a slight change in the Y expression with an additional factor of (3 —p) . 
These equations are solved same way as outlined in §4.1.1. We find that the solutions have the 
same dependence on z^^, & t-y as in equations llOlfflO^ And the dependence on r]i and rja are: 

_(p+18) 18-p 6-p 7p-18 36-17p p+18 

7i oc r]a ^® rij^^'^ , B (X {rjirja) 12 , T oc {rjirja) , and t (x r]^ rja ■ The distance of 7-ray 

7p— 18 . r) rj 

source (x T (x {rjir]a) ■ We constrain rji by requiring ~ ^^alt ~ 2^7??^ ^10 keV, which 
suggests that r/j<0.3, in accord with the numerical calculation. 

Numerically, we find 7i<20, B<2 x 10^ Gauss, r<103, and t>10-^ for GRB-\] (fig. [U]). The 
7-ray source radius Ry^2 x 10^^y~9 (l + y) 9 cm is in good agreement with R-yS5 x lO^^cm obtained 
by numerical calculations. Note that R^ for this case is a factor of a few higher than the Uc < t^i 
case discussed in §4.2.1. 

The numerical results for the allowed region of 5-D space, for z^c > t^i, are also shown in 
Figure [TTJ The solutions corresponding to Vi > Vc are those at larger radius - the right hand side 
of each solution contour, or bubble, shown in the figure. These solutions have smaller 7^, larger T, 
smaller B &i t, and a little bit smaller Y. Since these solutions too seem viable for shock models, 
we explore below the x-ray and optical flux accompanying 7-ray emission. 

4-2. 2a. Prompt X-ray and optical flux 

The analytical expression for the x-ray flux is almost identical to that found for the fc < 
case (eq. I106P - the only difference is in the dependence on r]a and r/j. We find that - 

(p-2)(p+18) 

fx OC (rjiTJa) 36 , (113) 

37 1 

and the lower limit on the x-ray flux is ~ BOY^se (l-|-y)^T8 mJy; this is in agreement with numerical 
results shown in fig. (jl2p . 

The R-band flux, dominated by synchrotron emission, is 



2eV\^ fiya\ 2 /i/A 2 



^-l^j [fj l^j *^"> 

which is, in terms of the observed quantities, given by 

o — i 5 1 2 2 2 _i 2 2 

~ 3 X 10^ u^,^f^H^tiY- HI + Y)sd^l{l + z)-3 Alp ir]aV^^— mJy (115) 

with the value of //j^l700 mJy for GRB-\\. This is higher than the numerical results that give 
/ij<200 due to the sensitivity of fR on r]a'r]i^2. If we assume that ta ~ t-y, as has been done 

numerically, we find (x z^^ ^/^t^. Numerically we find that fji increases with i/^, and t^, 
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and is most sensitive to t^. The increase of Jr with is again due to the high sensitivity on rjaTji 



which numericahy we find is oc i'- 



-1/2 



We note that the optical fiux for this case (fj < Vc) is larger than when Uc < i^i ~ 10^- lO'^mJy 
- since is larger for these solutions. For ta/t-y = 10~^ is reduced by a factor of ~ 50. This 
brings the highest flux for the ~ 1 s solutions down to about 140 mJy - still a bit bright, but 
there are many other solutions for this case that have flux smaller than about 10 mJy, and probably 
in accord with observations and upper limits. Thus, bright optical flux is a generic prediction of 
the SSC model for 7-ray generation whether the low energy spectral index is positive or negative; 
the optical is particularly brigh t for a > 0. This is a problem for the SSC model if 10^^ cm, 

as found in iKumar et al.l (j2007l ) , as the brightest optical flux is produced at the larger R^ . One of 
the, possibly only, ways to avoid bright, prompt, optical emission is if ta ^ t^. 



5. GeV photon signal for synchrotron and SSC solutions 



Detection of GeV photons by the Gamma Ray Large Area Space Telescope (GLAST) (jMcEnerv &: GLAST Mis 



20061 ) will be a useful piece of evidence to use to determine if GRBs are produced by synchrotron 
or SSC emission mechanisms. The IC scattering of 7-ray photons produced by synchrotron will 
peak above the GLAST band, ^ 100 GeV, while the second IC scattered SSC photons will peak 
at ~ 1 GeV. 

For the synchrotron cases a = 1/3 and —{p — l)/2, max(7j,7c) ~ 10^ and r<10^^, giving 
the peak of IC scattered flux i/fi, at z^g ~ min[zv^ max(7j, 7c)^, max(7j, 7c)rmeC^/(l + z)] 10^ 
GeV and flux ufi, ~ t/^ug ~ 10^^ erg s~^ cm^^. We need to ensure that such high energy 
photons can escape the source region and are not converted to electron-positron pairs. This effect 
is incorporated in our numerical calculations and is discussed below. Morever, photons of energy 
larger than ~ ITeV are converted to by collision with infrared photons and therefore we would 
not see ^ITeV photons from GRBs at cosmological distances. 

The SSC a > solutions have 7^ ~ 7c^lO^ and t>10~^ giving the second scattering peak of 
i/Q ~ 4z/^7? ~ 1 GeV and the flux z//jy ~ 10"'' erg s~^ cm~^. The SSC a < solutions have very 
similar z^g and ufi^. The SSC signals are well above the GLAST threshold so we expect to see a 
GeV signal from GRBs produced by the SSC process. For synchrotron solutions, however, the IC 
flux might be below the GLAST threshold. The spectral shape will also be different - is well 
above the GLAST band for the synchrotron case, while SSC should peak at the low end of the 
GLAST band. 

Using analytical results for the synchrotron a = 1/3 case (Equations H5ti48p , we find that the 
IC scattered signal peaks at 



ly'c ~ -fiTmeC^/il + z) ~ 2 X 10^4 4 i|y- 1 (1 + y)td|^^(l + z)-t 



P 



8 



GeV; (116) 
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This is due to the Klein-Nishina reduction to cross-section above electron rest frame photon energy 
of mgC^. 

The frequency above which the high energy spectrum is attenuated due to 77 — > e~e^ within 
the GRB source is i'± ~ {TrrieC^)'^ /uT-r^i, where i^r^i is the frequency of the synchrotron photon 
at which the optical depth to pair production with iy± photons is 1. The optical depth to pair 
production is r^-y ~ aTn'^^iR'y/T, where n'^, ~ Liso{i^)/4:'7TTcR'^ is the comoving number density 
of photons between v' and 2^', FLy/T is the comoving shell width, and Liso{i') is the observed 
isotropic luminosity per unit frequency. To find J^r^i, we set r-y^ ~ 1 and solve the equation using 
the observed 7-ray spectrum; i'± is calculated from T & Vrr^i- In terms of the observable parameters 
for the synchrotron a = 1/3 case, we find 

3 

' GeV (117) 

where we have assumed that the synchrotron GRB spectrum is Liso^v) oc z^^^ for v > v^. Since 
^± < (calculated above), and the spectrum falls off very steeply above i'±, the IC spectrum will 
peak at i/j- for many of the synchrotron solutions. 

The flux at Vq, with appropriate Klein-Nishina cross section, is 

ergs~^cm~^. (118) 

This flux is probably just at the GLAST threshold for detection. If i'± < Uq, the uf^, spectrum peaks 
at Z/-I-, and the flux at this frequency will be smaller than that in equation (jllSp : the attenuation 
of >TeV photons as they propagate through the inter-galactic medium would further reduce the 
observed flux. The results are very similar for the a = —{p — l)/2 case, since the a = 1/3 case is a 
subset of the a = —{p— l)/2 solutions with 7^ ~ 7c. 

For the SSC case, using the a > analytical results in Equations [20] - E31 the 2'^'^ IC peak is 

3 1_ 1 2 2 _i 2 — i P+18 

^}^^3ulfy'H~^tlYHl + Y)-9d^l{l + zy^A^;ria GeV (119) 
and the flux at this peak is 

[uUf ~ 7 X 10-' u^J^YA^^^ erg s'^cm'^. (120) 

I'G is smaller than i'± for the majority of the SSC solution space, so the 2nd inverse Compton 
scattering spectrum will indeed peak at i^g- For the SSC a < case, the expressions are very 
similar - the only difference is the dependence on t/q, rji, and p: 

3 _J_ 12, 2 — i 2 —1 P+18 18-p 

~ 34/^ (1 + yy^d^lil + z)-^A,;7]a ri, i« GeV (121) 

and the flux at this peak is 

[i^f.r ~ 7 X 10-' ly^JyYrjf-PA-^^ ergs-icm-2. (122) 



~ 4 X 10'uif^H-HaY-l{l + Y)dl(l + z)5 
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synchrotron: a = -(p-1)/2 synchrotron: a = 1/3 




10"^ 10"^ 1 10 10"^ 10"^ 10"^ 1 10 

[GeV] [GeV] 



Fig. 13. — Numerical results for the IC scattering of prompt 7-ray photons in the source. The four panels 
show the peak of IC spectrum (va) and the flux (vfu) at the peak when the 7-ray emission is produced via 
the synchrotron process (top two panels) and via the SSC process (bottom two panels) for a > & a < 
cases; for all of these cases we took the underlying 7-ray spectrum with i'^ — lOOkeV & = O.lmJy. The IC 
signal for many of the synchrotron solutions are affected by photon-photon pair production within the GRB 
source, and that is included in these numerical calculations; however the conversion of >lTeV photons to 
due to collision with infrared photons in the inter-galactic medium is not shown in the figure (no photon 
above 1 TeV can reach us from a GRB source at cosmological distances). The effect of pair production is 
very small in the SSC cases. The Klein-Nishina cross-section has been used in these calculations, and it 
significantly affects the IC peak flux when 7-rays are generated via the synchrotron process. 

We show numerical results for and vfu for 4 cases in Figure [13] - the two synchrotron cases 
of a = —{p — l)/2 and 1/3 and the SSC cases of a < and > 0. These numerical results are in 
rough agreement with our analytical estimates. 

In summary, IC scattering of prompt 7-ray photons, when the GRB emission is produced via 
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the synchrotron process, gives rise to a spectrum that peaks at ~ TeV, the flux is of order 10~ erg 
cm~2 (fig. 13), the spectrum below the peak is between i/^^^ or i/"^^^^)/^, and the spectrum 
above the peak is expected to be sharply cutoff due to pair production. On the other hand, if 
the GRB emission is produced via the SSC process then the spectrum of 2"*^ IC scattered photons 
should peak at ~ 1 GeV, with a flux of ~ 10"^ erg cm"^ s~^ (fig. 13), and the spectrum above the 
peak should be f,^ oc This signal should be detected by GLAST. If GRBs are not detected 

in the GLAST band, that would suggest that GRB prompt emission is not generated via the SSC 
process. We note that the GLAST band flux cannot be reduced in the case of repeated acceleration 
of electrons during a single GRB pulse. 



6. Comparison with prior work on 7-ray generation mechanism 

W e provide a brief comparison with published work on 7-ray generation in the int ernal shock 



moc 



de l dRees fc Meszarosll 19941 : IPaczvnski fc Xulll994l : iPapathanassiou &: Meszaroslll996l : lMeszaros &: ReesI 



19971 : ISari &: Piranlll997bl ) and the formalism/results of this paper. There is a fairly extensive lit- 
erature on this topic and this is not the place to provide a general review. What we wish to do is 
to describe the main difference between previous approaches and the work presented here. 

We should note that 7-ray generation in the external shock model has also been looked at by 
a number of people eg. Rees & Meszaros (1992), Meszaros & Rees (1993), Piran et al. (1993), 
Dermer et al. (1999), Dermer & Mitman (1999), McMahon et al. (2004), Ramirez-Ruiz & Granot 
(2006); the issue of variability in external shocks is discussed in Sari & Piran (1997a), Dermer &: 
Mitman (1999), Nakar & Granot (2006). We don't have anything particularly enlightening to say 
regarding the external shock model that has not already been mentioned by one of these authors; 
the general problem with shocks is discussed in §3 & §4. 

The main difference between previous works and our approach is that previous works considered 
the forward problem i.e., starting with a parameterization of the properties of colliding shells and 
resulting shocks the emergent radiation field was calculated, whereas our approach is to start 
with the minimum number of physical parameters needed (five) to calculate the observed fiux 
and spectrum - at one instant in time or for one pulse in a multi-pulse GRB lightcurve - and 
determine these using the observed data. For synchrotron & IC radiations the parameters needed 
are 7^, F, B, N, and r, which are determined by the observational data u^, t^, and a for a pulse 
in GRB LC. The 5 parameters in turn are used to provide constraint on the nature of GRB source. 
The old forward approach is wedded to a particular model - either internal or external shock - 
whereas the method used in this work is relatively model independent. 

The parameters of the internal shock model can be mapped into the five parameters (7^, F, 
B, N, t) in a straightforward manner - this in fact is done implicitly in all the forward approach 
papers in order to calculate the emergent radiation. The converse of this is not true, however, since 
the internal shock model has more than 5 independent parameters (appendix A describes how to 
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go from the five parameters to providing a limit on some of the internal shock parameters such as 
the initial LFs of colliding shells and their comoving densities). 



Papathanassiou &: Meszarod (|l996l ) and lSari &: PiranI (jl997bl ) carried out a fairly detailed anal- 



ysis of prompt emission in the internal shock model. These authors addressed a set of questions 
such as the ability of synchrotron/SSC in internal shocks to produce a spectral peak near 100 keV, 
the observed flux in the 7- ray ba nd, and short time scale variability. iPapathanassiou &: Meszaros 



(|l996l ) and ISari &: PiranI (jl997bl ) realized that the cooling t ime scale for electron s (compared to 
the dynamical time) for internal shocks is short a nd although ISari &: PiranI (Il997bl) d on't explicitly 
say, their work applies to GRBs with a = —1/2. IPapathanassiou &: MeszarosI (|l996l ) look at com- 
posite synchrotron/SSC spectra, however there are many free parameters and little comparison to 
observed properties of GRBs. Had these authors investigated the self-consistency of synchrotron 
internal-shock model for the case of a = 1/3, they would have discovered the problem reported in 
this paper using their forward modeling approach. 

Ghisellini et al. (2000) did in fact worry about synchrotron solutions when a > and concluded 
that it is impossible to account for it in shock based models (this is paraphrasing their actual 
wordings). They pointed out that electron c ooling time, during prompt e r nissio n, is much smaller 
than th e dynamical time [as was reported in Papathanassiou Sz Meszaros (jl996l ) and Sari Piran 



( 1997bl )] and therefore the GRB spectrum below the peak should be always if the radiation 

is produced via synchrotron process. Ghisellini et al. (2000) considered the possibility that a lower 
strength magnetic field would avoid the excessive cooling of electrons, and rejected it based on 
the argument that IC emission would dominate in this case i.e., y ^ 1, and that IC spectrum 
too would be falling of as or steeper due to IC cooling. We find that a smaller magnetic 

field can avoid excessive cooling, so that a = 1/3, and at the same time Compton y ~ 1. The 
reason for these different conclusions is that we don't impose any restriction on the source distance 
that forward calculations based on internal shocks do. The most serious problem with synchrotron 
a = 1/3 case is that i?-y > Rd unless no <C 1 cm~^ (see §3.2). 

Ghisellini et al. (2000) correctly pointed out that re-acceleration of electrons, in shock based 
models, would not work because it requires too much energy; one has a continuous stream of 
electrons crossing the shock front ~ and to accelerate all of these electrons to their original energy 
distribution, so that a = 1/3, while they are rapidly losing energy to radiation will indeed require 
much more energy than is available in the shock. The re-acceleration invoked in this work is not 
shock based. It in fact requires abandoning shock models and considering a scenario where particles 
are NOT being added to the "system" - the source for 7-ray photons - continuously (as in shock 
based models) but where the source has a fixed number of particles that are being continuously 
accelerated; there is no excessive energy problem in this scenario. 
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7. Discussion 

In this paper we have investigated the generation of 7-rays in gamma-ray bursts via synchrotron 
or synchrotron-self-inverse-Compton (SSC) emissions in a relativistic outflow. The SSC radiation 
from a relativistic source can be fully described by a set of 5 parameters (7^, T, B, N, r); see table 
1 for definition of symbols used in the paper. For each possible low energy spectral index, we have 
analytically and numerically determined the region of the 5-D parameter space that is consistent 
with a set of GRB observations - i'^, f^, & t^. For these allowed regions - or solution sub-space - 
we calculate the x-ray and optical fluxes that should be seen concurrent with the 7-ray radiation 
to further narrow down the properties of 7-ray sources. The set of five parameters also allows us 
to determine the distance of the source from the center of explosion (Ry). 

We find that if 7-ray emission were to be produced via the synchrotron process, the required 
set of parameters and burst radius have extreme values that are not internally consistent and are 
in conflict with afterglow data. In particular, when the low energy spectrum is f^, cx u^^^ or u 2~ , 
the Lorentz factor of the source is required to be larger than ~ 10^, in disagreement with afterglow 
modeling, and the source distance (-R-y) is larger than the deceleration radius even when the density 
of the medium is as small as ~ 0.1 cm~'^. The requirement on the magnetic field strength is also 
very stringent; the comoving field strength is required to lie in a very narrow range of about ~ 10-30 
Gauss in order to explain the radiation for a typical burst as synchrotron emission. Allowing for the 
possibility of multiple electron acceleration episodes, i.e. ta <^ t^, alleviates the problem of large 
&: F; in this scenario the synchrotron process could account for the prompt 7-ray radiation for 
GRBs (see §3.1 & 3.2) although is still larger than the source distance determined for a subset 
of bursts detected by Swift (Kumar et al. 2007). 

The reason that synchrotron solutions with a = —{p—l)/2 &: 1/3 require large R^ & F is easy 
to understand. The number of electrons needed to produce the observed flux via the synchrotron 
process is ~ 10^^/-y/(SF). And in order to keep the Compton-Y parameter {Y ~ r^i^c) less than 
~ 10 — otherwise most of the energy will come out in IC-scattered photons at S> 1 GeV — the 
source must have small r or large R^. Moreover, since ~ -R^/(4cF^), large R^ solutions also have 
large F for a given GRB pulse width of t^. The reason that ta <C offers a way out of this problem 
is also easy to understand. Frequent re-acceleration of charge particles makes it possible to have 
larger magnetic field while keeping fc^lOO keV. This decreases the number of particles required to 
produce the observed flux (/^), and that in turn makes it possible to have a smaller R^. 

For Jy (X the allowed region of the 5-D parameter space for the synchrotron process 

is quite reasonable. However, interpreting these solutions in terms of the internal shock model 
requires the ratio of LFs of the two colliding shells to be rather large (>10 — 20) when the ratio of 
magnetic energy to electron energy is ^10, i.e. if we want the energy fraction in electrons to be not 
too small - otherwise 7-ray production would be very inefficient (see §3.3). 

The SSC process provides viable solutions for the prompt emission of a large fraction of GRBs. 
We have considered almost all different possibilities of the low energy spectrum for GRBs: fl'^ oc 1/°' 
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below the peak of I'fl'^ with — l^a < 1. The solution space (a hypersurface in the 5-D parameter 
space) is quite large for a < and 0.5<a < 1. However, there are no SSC solutions when a ~ 1/3; 
the reason is that the synchrotron characteristic and cooling frequencies should be equal {vi ~ Vc) 
in order that a = 1/3, and in that case the synchrotron-self-absorption frequency is shown to be 
roughly equal to Ui as well (see §4.1.1), and therefore the low energy spectral index is ~ 1 and NOT 
~ 1/3 as desired. 

The SSC solution space for a < and 0.5^a < 1 has source LF of order 100, the minimum 
electron energy <10'^meC^ (characteristic of mildly relativistic shocks), and 10^^<i?-y<10^^cm is 
smaller than the deceleration radius. These solutions are accompanied by bright optical synchrotron 
flux of ~ 10 mJy (14*^ mag) to several hundred mJy for bursts at z ~ 2 - brightest for bursts 
with a ~ 1 and those bursts with pulse duration on the order of >ls. Moreover, the optical flux 
is correlated with - the flux is larger for larger - and for i?^>2 x 10^^ (cf. Kumar et al. 
2007), the optical flux is ^lOOmJy (< 12*^^ mag). Bright optical flux contemporaneous with 7-rays 
is a prediction of the SSC model that is in conflict with prompt optical follow up observations of a 
large number of bursts detected by Swift (Roming et al. 2006) |j We note that the large optical flux 
accompanying 7-rays can be reduced only if each radiating electron is accelerated numerous times 
(>100) in time period of order the duration of a pulse in the GRB lightcurve, i.e. ta^t^/W'^- GRB 
models based on converting kinetic energy to radiation via shocks have ta ~ t^, where particles 
are accelerated at the shock-front and not down-stream, but continuous particle acceleration might 
work for some alternate scenarios such as magnetic field reconnection/dissipatio n. Ma g netic outflow 



Thompson ( 


1994). 


Katz ( 


1997) 



& Konigl (2001), Spruit et al. (2001), Lyutikov & Blandford (2003), Thompson (2006). However, 
the model has not been developed to the extend where it can be tested with GRB observations. 

The data from the upcoming high-energy mission GLAST should be able to settle the question 
whether GRBs are produced via synchrotron or the SSC process (see §5); see Gupta & Zhang 
(2007), Granot et al. 2007 (and references therein) for recent work on how GLAST would help our 
understanding of GRBs. 

To summarize our main conclusions, we find that the synchrotron process has serious difficulty 
accounting for the prompt emission in GRBs. The SSC offers reasonable solutions for all GRBs 
except those with spectral index of 1/3 below the peak. SSC solutions predict very bright optical 
emission (> 10m Jy or 14-mag for 2: ~ 2) accompanying 7-ray lightcurves which is in conflict with a 
number of well observed bursts. A possible solution to this problem might be to drop the assumption 



^The simultaneous optical and 7-ray observations for a few bursts show bright optical flares. For example, GRB 
041219a was observed in optical and IR simultaneously with 7-rays jVestrand et ahlboosl : iBlake et allboosl ): the 
optical flux peaked at 13.7 magnitude at approximately t he same tim e as th e first main pulse of the GRB lightcurve 
(as expected for the SSC model), and IR measurements iBlake et al. show evidence of rapid variability with 

the last spikes in the GRB light curve. GRB 990123, with a positive a, had a peak optical flux during the burst of 



IJy jAkerlof et al.lll999l ). however the optical lightcurve was not correlated with the 7-ray lightcurve. 
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that ta ^ try. The assumption of one shot acceleration of electrons i.e., ta ~ t^^, is motivated by 
shock based physics for GRBs and it may have to be replaced with an alternate scenario in which 
all electrons that radiate in the 7-ray band are accelerated continuously throughout the duration 
of a 7-ray pulse. In that scenario there are viable synchrotron solutions when a ~ 1/3 - a case 
that otherwise cannot be explained by the SSC mechanism. 
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A. Determination of the relative Lorentz factor of colliding shells 



We can relate solutions (7j, F, B, N, r) we find for a set of observables to parameters for any 
model for GRBs. In this appendix we relate our solutions to the parameters for the internal shock 
model where 7-rays are produced by collision of two shells, as shown in Figure [T^ Internal shocks 
are produced by a collision between fast ejecta catching up with slower ejecta (IRees &: Meszaros 
19941 : iPaczynski &: Xulll994l ) which in the discrete version is modeled as the collision between two 
homogeneous shells moving with LFs Fi and F2 (Fi is LF of the faster, inner shell). [The external 
shock can be thought of as a special case of internal shocks which results from collision between a 
stationary circum-stellar medium (F2 = 1) and the ejecta from the burst.] 



shock front 2 — 1 




stiock front 1 



Fig. 14. — Pictorial representation of two colliding shells — internal shock model. 

This appendix is devoted to describing the method we use to place a lower limit to the relative 
LF of the colliding shells from T^h, the LF of the shock front in the 7-ray producing shell, which 
is directly related to 7^ (one of the 5 parameters) via equation pT|) . and by varying the ratio of 
densities in the inner and outer shells (ni/n2) subject to the condition that the efficiency for 7-ray 
production is not less than ~10%; as a by product we find the comoving density ratio for the shells. 



The LF of unshocked inner shell, which is moving faster and lies a bit closer to the center of 
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explosion than the outer shell, can be determined from the addition of LFs Tgh and F - 

Ti = r,hr{i + vshv), (Ai) 

provided that the 7-ray emission we observe is produced in the inner shell when it is shock heated 
{v &: Vsh Sire speeds corresponding to LFs F - the LF of the 7-ray source - & Tsh)- On the other 
hand if it is the outer shell that produces the observed 7-ray emission then its LF factor, before it 
was shocked, is given by 

F2 = F,^F(1- w). (A2) 

The LF of the shell that does not contribute significantly to the observed 7-ray emission cannot 
be determined uniquely. The problem is that in the absence of any observed emission that can be 
identified with this shell we cannot say anything directly about the LF of the shock front moving 
into this shell. However, we can still constrain its LF by requiring that the comoving density of this 
shell is such that the efficiency for thermal energy produced in the collision is no less than ~10% 
(most GRBs for which we have good afterglow data and we can determine the kinetic energy of 
the ejecta show that the efficiency for converting kinetic energy to 7-rays is about 10% or more). 

The efficiency of 7-ray production is used to constrain ni/n2, and that in turn provides a lower 
limit for the relative LF Frei ■ The fraction of kinetic energy converted to thermal energy when two 
shells of mass mi and m2 collide with a relative LF of Frei is (Kumar, 1999; Piran, 1999) 



^ _^ 2{m2/mi){Trei - 1) 
[1 + m2/mi]2 



(A3) 



For equal mass shells a Trei = 1-5 collision converts 10% of the kinetic energy of shells to thermal 
energy, and for Frei > 5 the conversion efficiency is more than 40%. We need, however, to figure 
out the fraction of energy produced in a collision that is radiated in the typical 7-ray observing 
band of 15-400 keV during the time interval tj, in two shells that are not of equal mass. 

The ratio of the mass of the two shells is the mass swept up by the two shocks in the time 
equal to the shock transit time for the shell that produces the observed 7-ray emission. This is 
given by 

(A4) 



mi ni{4rs-,+3)vss-, 



m2 n2{4Ts,2 +3)vss2 

where rii and n2 are the densities of shell 1 & 2, F^j and F^j are the shock LFs in the frame of each 
unshocked shell (if the GRB emission is predominantly from shell 1, then F^j = F^/j), and Vgg^ and 
Vss2 are the speeds of the shock fronts with respect to each unshocked shell. 

The shock front speeds are determined from the following cubic equation obtained from the 
continuity of energy, momentum and particle number fluxes across the shock front (Landau &i 
Lifshitz, 1980): 

rL + fi - r V-rL - f 1 - A V- - f 1 - -1 =o. (As) 



where Tgu & Tgs are the LF of the shock front wrt to the unshocked and shocked fluid respectively, 
and a depends on the equation of state of the shocked gas and is approximated by 

4ro + i 

i + i 

which provides a smooth interpolation between the sub-relativistic value for = 5/2 and the highly 
relativistic value of 4. 

The shock LF wrt the unshocked fluid in the shell not dominating the GRB emission (Fsj, if 
we assume that the GRB is being produced in shell 1) can be determined from the condition of 
pressure balance across the surface of discontinuity separating the shocked fluids in the two shells- 

ni (4r,, + 3) (F,, - 1) = n2 (4F,, + 3) (F,, - 1) . (A7) 



For a given density ratio ni/n2 and F^^ we can solve the above equation to determine F^j, 
which in turn is used to determine shock front speed wrt the unshocked fluid for the outer/inner 
shell using 

and the cubic equation IA5l for F^gj. These pieces together give us mi/m2 and Trei, which are used 
to calculate - the fraction of the kinetic energy of the two shells converted to thermal energy in 
shell collision (eq. A3). 

The ratio of the internal energy of the shocked gas in these two shells is 

Si ^ ni(4F,, +3)(F,, 
E2 n2{ATs,+2,){Ts^-l)vss2 

The fraction of the total internal energy of the shocked gas in shell '1' is then /i = -g^/ |^-^ + ij . 
We assume that the majority of this energy is indeed being radiated in the GRB band. For shell 
'2', /2 = 1 — /i, and we find the fraction of the total radiation contributing to the GRB band, 
fcRB] fcRB is found from the ratio Ui^U^ (t^h) / 1^12 (1^12 )> where Uiji^i^ = (721/7*2)^ (assuming 
that the magnetic field is equal in both shells). The radiation efficiency in the GRB band for the 
shell collision is then Jr (/i + f2fGRB)- 

When considering synchrotron radiation as the primary source of emission in the 7-ray band we 
need to take into account the energy fraction that is lost to very high energy photons (z/ ^ 200 keV) 
produced via the inverse-Compton process. The fraction of energy radiated via the synchrotron 
emission is 1/(1 -|-y), therefore, the total efficiency for energy production in shell collisions must be 
larger than the desired 10% by a factor of 1 + Y - this effectively restricts solutions to y^lO. y^O.l 
if inverse-Compton emission is the main source for the observed 7-ray emission and also Y must not 
be greater than ~ 10^ otherwise most of the radiative energy will be in the 2nd inverse-Compton 
photons of much higher energy. In the same sense, we need to make sure the ratio of magnetic 
energy to that in electrons is <1, in order for the electrons to radiate efficiently. 



(A9) 
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Fig. 15. — Left panel: Dependence of the ratio of tlie two sliell masses, toi/to2, on tlie ratio of densities, 
nxjn^. The colored lines show the relation for different relative Lorentz factor, r^ei, as shown in the legends. 
The dotted and dashed lines show the applicable dependence on the density ratio. For high V^ei and 

1 /2 

nxjui > 0.01, mi/m2 oc (ni/n2) , however the deviation from this relation for small Trei and ni/n2 is 
small - about a factor of 2. Right panel: Dependence of the ratio of synchrotron injection frequencies (vi) 
in the two shells on ni/n2. For high Trei and ni/n2 > 0.01, i'isl^i,2 oc n2/ni. At low Trei and small ni/n2, 
the ratio of injection frequencies in the two shells significantly deviates from this, and for mildly relativistic 
Trei, can be better approximated by Vi^i/vi^^ oc (n2/ni)^, shown by the the dashed line. 

For highly relativistic forward and reverse shocks Vssx ~ Vss2 ~ 1/3, ~ F^^ (ni/n2)^/^, 
Fre/ ~ FsjFs2, mi/m2 ~ (?^i/f^2)"^^^, and the ratio of the characteristic synchrotron frequency 
in shell 1 to shell 2 (assuming the same magnetic field strength) is ~ n^lnx. We note that the 
assumption of highly relativistic shocks is not valid for many solutions we find for the prompt 
7-ray emissions, and that all the results reported in paper are obtained by numerically solving the 
appropriate equations. The numerically solved relationships of the ratio of the masses and injection 
frequencies in the two shells are shown in Figure [T5l and compared to the analytical estimates for 
highly relativistic shocks. 

In summary, for a given Fs/i(7i) we vary nxjn^ and determine the mass ratio and the relative 
LF of the shell collision (F^e/) so that the gamma-ray production efficiency is above a certain desired 
value (10%). All of the numerical results we show for V^ei were calculated using these steps. Using 
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our upper limit on we also calculate the expected emission from the shell 2 in the x-ray and 
optical bands, and include this in our analysis. 



B. Jitter radiation process and GRBs 



A radiation mechanism, called jitter radiation, has been suggested by iMedvedevI (|200Cl ) as 
the process for 7-ray emission for those cases where the low energy spectrum rises more steeply 
than v^l"^ expected of the synchrotron radiation (such spectra are said to lie "above the line of 
death" because a non-self-absorbed synchrotron spectrum cannot have this steep of a rise). The 
jitter radiation is produced when the coherence length scale for magnetic field is short and electron 
trajectory is perturbed before it has traveled a distance of a Larmor radius. This is an attractive 
idea for explaining a class of GRBs lying above the line of death, and we explore its applicability 
to GRBs in this appendix. 

The peak jitter frequency in lab frame (as seen by an observer at rest in the host galaxy) is: 

Vj = X^'Kq^V shUel (me7e)7pr, (Bl) 

where q is electron charge, F is the bulk Lorentz factor (LF) of the source, 7^ = min{7j,7c}, 7c 
is the thermal LF of electrons that cool on a dynamical time, 7^ is the minimum thermal LF of 
electrons behind the shock front (note that 7j = ee{mp/me)Tship — 2)/(p — 1); where Tgh is the LF 
of the shock front wrt the unshocked shell), rig is the comoving electron density in the unshocked 
shell, and 7e ~ 3 - 4 is the initial effective thermal Lorentz factor of the streaming electrons. 

In order to get the gamma-ray burst spectrum below the peak to be proportional to ~ j^^, we 
want Uj ~ 10^ eV (the peak of the GRB spectrum is of order 100 keV). Therefore, from the above 
equation we find the following condition on comoving electron density in the unshocked shell: 

4 X 10^^ 

^ pay i%/'^sh) (B2) 



or the optical depth of the source to Thomson scattering is: 

T « An,T,hR^aT/T « \ ^ (B3) 

r7p 



where i?^ is the distance of the source from the center of the explosion, and ~ i?^/(4F^c) is the 
GRB variability time scale. 

For internal shock V is of order a few, and therefore, 7^ ~ 10^. In this case we find the optical 
depth of the source to be 

r-i^. (B4) 
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Or r ~ 1 for ~ 10 ^ sec, and F ~ 10^. The next step is to estimate the LF of cooling electrons 
(7c), which can be shown to be: 

67rmec(l + z) 
^''^ aTB%T{l + Yy 

where B is the magnetic field in the source co-moving frame, and Y is the Compton-Y-parameter; 
Y Tjc7i (for p^2). 

If 7c is not much less than ji then Y is very large (of order a million) , and most of the radiation 
energy will come out as SSC photons at > 10^ GeV. And moreover, electrons cool very rapidly 
resulting in the synchrotron cooling frequency (fc) to be much less than 100 keV, and therefore the 
spectrum below the peak will be and not cx u^. Even if we take 7c to be order unity, we still 

run into similar problems. 



We, therefore, do not include the jitter process in our analysis of GRB prompt radiation 
mechanism. 



